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Estuarine  Circulations  and  Mass  Distributions 


1.  Introduction 


These  lectures  concern  estuaries  or,  more  precisely,  my  special  recent 
interests  in-the  physical  aspects  of  estuaries.  Donald  Pritchard  at  The  Johns 
Hopkins  University  has  studied  these  systems  for  twenty-five  years  and  has 
attempted  (1967)  to  formulate  a workable  definition: 

*An  estuary  is  a semi-enclosed  coastal  body  of  water  which  has  a free 
connection  with  the  open  sea  and  within  which  sea  water  is  measurably  diluted 
with  fresh  water  derived  from  land  drainage.  * 

J 

I can  quarrel  with  this  definition  in  certain  respects.  Pritchard  apparently 
wanted  a definition  which  required  an  estuary  to  have  an  "estuarine  circulation" 
characterized  by  ' /ater  outflow  in  the  upper  layercomposed  of  a mixture  of  fresh 
water  and  salt  water  and  an  inflow  in  the  lower  layers  of  saltier  sea  water1  . The 
presence  of  this  kind  of  circulation  distinguishes  an  estuary  from  other  types  of 
embayments  and  is  a useful  characterization.  However,  as  Pritchard  himself 
states,  it  excludes  a body  of  water  such  as  the  Baltic  Sea  by  including  the  adjective 
"coastal".  The  Baltic  has  an  estuarine  circulation  and  I think  it  should  be  considered 
an  estuary.  A further  objection  is  to  the  last  four  words  in  the  definition  which  put 

an  undue  importance  on  the  source  of  the  fresh  water.  It  may  be  interesting,  for 
example,  to  consider  a theoretical  model  of  an  estuary  in  which  the  direct  addition 

of  rain  water  may  be  an  important  source.  Certainly  we  should  not  exclude  the 


1 In  the  typical  estuary  density  differences  due  to  salinity  are  considerably 
greater  than  those  due  to  temperature  (Pritchard,  1965). 
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fresh  water  addition  from  the  rain  falling  directly  on  a body  as  large  as  the 
Baltic  Sea. 


Pritchard's  definition  of  an  estuary  requires  an  appreciable  amount  of 
incoming  fresh  water.  However,  we  will  see  that  interesting  three-layer  circu- 
lations with  resemblances  to  two-layer  estuarine-type  circulations  can  occur  when 
the  fresh-water  supply  is  zero  or  very  small,  and  we  will  include  this  case.  We 
may  even  go  further  and  include  situations  in  which  evaporation  exceeds  precipi- 
tation and  run-off.  In  such  cases,  for  example  the  Mediterranean  Sea  and  Laguna 
Madre  in  Texas  (Pritchard,  1965),  inverse  two-layer  circulations  exist  which 
also  resemble  the  classic  estuarine  circulation. 

Pritchard  goes  on  to  distinguish  several  types  of  estuaries  based  on  geo- 
logical considerations  of  the  origin  of  the  estuaries,  for  example  drowned  river 
valleys,  which  he  calls  "coastal  plain  estuaries"  and  coastal  indentures  gouged 
out  by  glaciers  which  he  calls  "fjords".  Certainly  fjords  differ  fundamentally  from 
such  coastal  plain  estuaries  as  the  Chesapeake  Bay  but  the  way  in  which  the 
estuaries  formed  is  not  directly  important  for  us. 

It  will  clarify  our  problem  if  we  attempt  to  classify  estuaries  in  a manner 
similar  to  that  of  Stommel  (1951),  based  in  part  on  the  relative  strength  of  the 
turbulence  that  causes  the  mixing  of  waters  of  different  salinity  (density).  We 
define  four  types : 


river  emptying  into  the  ocean 


the  absence  of  any  mixing  the  upper  layer  remains  fresh  and  the  mass  flux  across 
any  section  is  equal  to  the  river  discharge  R.  There  is  a zero  salt  flux  across 
every  section.  Since  the  interface  is  stationary,  if  the  lower  fluid  may  be  con- 
sidered frictionless  and  barotropic,  it  will  be  at  rest.  Among  other  things,  this 
means  that  the  horizontal  pressure  gradient  in  the  lower  layer  will  be  zero.  The 
salt  wedge  may  penetrate  up  the  river  50  km  or  more  (Pritchard,  1965). 

Of  course,  frictional  effects  and  mixing  are  never  entirely  absent.  For 
example  even  in  the  absence  of  wind  and  tide  the  flow  of  the  river  water  over  the 
salt  wedge  will  tend  to  produce  unstable  waves  at  the  interface  which  will  break 
and  cause  mixing.  This  has  been  studied  by  Keulegan  (1949)  in  a laboratory  flow 
of  a light  fluid  over  heavier  salty  fluid.  As  pictured  in  Fig.  1.2,  waves  formed 
and  eddies  were  ejected  from  the  crests  causing  some  of  the  lower  salt  water  to 
be  mixed  with  the  upper  fluid.  In  the  river  this  will  cause  the  upper  fluid  to  become 
somewhat  salty  as  it  flows  toward  the  sea.  The  salty  water  comes  from  below  so 
that  there  must  be  a (small)  flow  in  the  lower  layer  from  the  sea  up  the  estuary  to 
replenish  the  lost  salt.  In  this  kind  of  mixing  none  of  the  upper  fluid  is  mixed 
into  the  lower  fluid  and  the  lower  water  remains  of  the  density  of  water  in  the 
sea.  Keulegan  found  that  mixing  began  when 

e = /Q  = 0.178 


where  vis  the  kinematic  viscosity,  g is  gravity,  Ap  is  the  density  difference, 
pr  is  the  density  of  fresh  water  and  u is  the  velocity  of  the  upper  fluid.  For 
typical  values  in  a river,  this  corresponds  to  velocities  of  a few  centimeters 
per  second  so  that  instabilities  of  this  kind  are  very  likely.  The  definition  of  a 
Type  1 estuary  is  restricted  to  situations  in  which  the  (weak)  mixing  is  primarily 
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due  to  the  wave  breaking  at  the  interface.  An  example  of  a Type  1 estuary  is  the 
mouth  of  the  Mississippi  River. 

Type  2 Estuary. 

Two  other  effects  are  important  in  producing  mixing  in  estuaries.  The  first 
is  the  tide  which  produces  motions  and  breaking  internal  waves  in  the  fluid  and 
these  in  turn  cause  turbulence  and  mixing.  The  second  is  the  stress  of  the  wind. 
This  causes  velocities  in  the  water  and  the  resulting  vertical  shears  are  a source 
for  energy  of  turbulence.  The  result  is  increased  mixing  compared  to  the  Type  1 
estuary  and  a density  profile  resembling  that  in  Fig.  1.  3a  with  a fairly  continuous 
density  distribution  with  height.  Here  the  density  increases  gradually  with  depth 
although  there  is  frequently  an  identifiable  halocline  across  which  the  density 
increase  with  depth  is  more  marked.  Examples  are  the  Chesapeake  Bay  and  the 
mouth  of  the  James  River  which  empties  into  the  bay. 

We  include  in  Type  2 the  fjord-type  estuary  which  is  typically  so  deep  that 
the  lower  layers  are  relatively  inert  with  little  mixing  and  a density  equal  to  that 
of  sea  water.  In  this  case  a relatively  sharp  halocline  often  exists  as  shown 
in  Fig.  1.3b. 

The  mixing  is  typically  much  larger  in  Type  2 than  in  Type  1.  In  the 
Type  1 the  flux  of  salt  water  into  the  upper  layer  is  less  than  or  at  most  of  the 
order  of  the  river  flow.  In  Type  2 the  upward  flux  of  salty  water  is  much  larger 
so  that  the  total  flow  of  brackish  water  seaward  in  the  upper  layers  will  be  several 
times  the  river  flow.  In  the  James  River,  for  example,  the  total  discharge  is 
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Type  3 Estuary. 

If  the  source  of  mixing  is  strong  enough  the  estuary  may  be  sectionally 
homogeneous,  i.  e. , thoroughly  mixed  from  top  to  bottom  and  side  to  side.  Al- 
though the  vertical  and  lateral  gradients  of  salinity  vanish,  a horizontal,  seaward 
variation  of  salinity  may  exist.  In  the  idealized  case  there  is  no  estuarine  type 
circulation  in  the  estuary.  The  observed  longitudinal  gradient  of  salinity  requires 
a transfer  of  salt  up  the  estuary  and  this  must  be  accomplished  by  turbulent 
diffusion. 

It  is  possible  that  strictly  sectionally  homogeneous  estuaries  do  not  really 
exist.  For  example,  the  estuary  of  the  Mersey  River  in  England  was  thought  to 
be  vertically  homogeneous  until  Bowden,  Fairbairn  and  Hughes  (1959)  reported  a 
slight  vertical  salinity  difference.  This  raises  the  basic  question  as  to  the  relative 
importance  of  the  advective  salt  flux  SS  set  up  by  the  estuarine  circulation  and  the 
turbulent  longitudinal  flux  u'S'.  In  these  expressions  u is  velocity  and  S is  salinity. 

If  u = 0,  the  estuarine  circulation  is  absent  but,  because  S is  large  compared  to  S', 
a very  weak  circulation  can  transport  a great  deal  of  salt. 

Type  4 Estuary. 

In  this  estuary  the  fresh-water  influx  from  river  discharge  is  negligibly 
small.  If  the  water  in  the  sea  is  of  homogeneous  density,  it  will  fill  the  estuary 
and,  since  no  density  variations  exist,  there  will  be  no  gravitational  circulations. 

The  Type  4 estuary  can  have  motions,  however,  if  there  is  a vertical  gradient 
of  density  in  the  outside  water  and  if  the  mixing  is  greater  in  the  estuary  than 
in  the  outside  water.  Then  we  will  have  the  situation  portrayed  in  Fig.  1. 4. 

The  constant  density  surfaces  slope  as  shown  and  this  leads  to  accelerations 
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which  produce  the  three-layer  circulation  pictured  in  the  figure.  This  possibility 
was  first  raised  by  Hachey  (1934).  A circulation  of  this  kind  was  inferred  from 
a study  of  acid  wastes  in  the  water  in  Baltimore  Harbor  (Stroup,  Pritchard,  and 
Carpenter,  1961).  The  Chesapeake  Bay  outside  the  harbor  has  the  required  density 
variation  with  height  caused  by  the  presence  of  lighter  water  in  the  upper  layers 
originating  in  the  Susquehanna  River. 

This  concludes  the  introduction  to  these  lectures.  In  subsequent  sections 
we  will  consider  some  aspects  of  estuaries  in  detail  based  on  recent  research. 


2.  Salt-Wedge  Estuary 

Let  us  consider  first  an  estuary  of  Type  1 pictured  in  Fig.  1. 1.  We  will 
confine  attention  to  an  idealized  situation  in  which  we  may  neglect  mixing  entirely. 
Then  the  lower  layer  is  at  rest  and  therefore  the  longitudinal  pressure  gradient 
along  the  estuary  is  zero  in  the  salt  water.  From  hydrostatics  (and  using  the 
Boussinesq  approximation)  the  pressure  in  the  lower  layer  is 

p/pr  = g(H-z)  + &>(H-h-z)  (2.1) 

where  b is  buoyancy,  defined  by  b = g(p-pf)/pf  and  is  the  buoyancy  of  the  salt 
water.  The  equation  of  the  free  surface  is  z = H,  and  h is  the  thickness  of  the 
upper  layer.  Thus  in  the  lower  layer 

gHx  =hx£b  (2.2) 

or 

gH  = hdb  + const  (2. 3) 


We  assume  irrotational  motion  in  the  upper  layer  so  that  the  Bernouilli  equation  is 

P us 

-L-+  + gz  = const 

Pf  2 


(2.4) 
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where  u is  the  longitudinal  velocity  which  will  be  nearly  uniform  across  a section. 

We  neglect  the  other  two  components  of  the  kinetic  energy.  They  will  be  very 

small  if  the  length  scale  along  the  estuary  is  large  compared  with  the  depth  and 

if  the  change  of  width  is  gradual.  Hydrostatics  leads  to 

a 

gH  + ^ = const  (2.  5) 


or 

h 1 + | y = E (2.  6) 

where  E is  a constant  proportional  to  the  sum  of  kinetic  and  potential  energies,  and 
F is  the  densimetric  Froude  number.  It  may  be  written 

^ = WV  to  (2,7) 


where  R is  the  river  discharge  and  W is  the  width.  Let  us  non-dimensionalize 
using  the  depth  of  the  fluid  ho  at  the  tip  of  the  salt  wedge  and  the  width  W0  of  the 
estuary  at  this  section.  We  have 

2h'  + = E'  (2.  8) 

(h-WT 


where 


W 


R3 

W^hgAb 


(2.9) 


and  E'  is  non-dimensional  constant  proportional  to  the  energy.  If  energy  is  con- 
served along  the  flow  E'  = 2 + F§.  Differentiating  (2.  8),  we  get 

= U^/d-F2)  (2.10) 

This  shows  that  the  depth  of  the  upper  level  decreases  with  distance  downstream 
(assuming  the  width  increases  downstream)  if  and  only  if  conditions  arc  super- 
critical, i.e. , F3  >1.  Eq.  (2.6)  shows  that  F3  then  increases  along  the  estuary. 
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Thus,  if  we  require  energy  conservation  and  that  h decreases  everywhere  toward 
the  mouth , Fq^  > 1,  i.  e.  conditions  must  be  supercritical  at  the  top  of  the  salt 
wedge  and  therefore  supercritical  everywhere.  In  fact,  it  is  physically  obvious 
that  h must  decrease  ultimately  as  we  proceed  seaward  so  that  conditions  must 
ultimately  be  supercritical  in  any  case.  It  is  possible,  however,  that  h could 
increase  near  the  tip  for  some  distance  so  that  we  may  consider  the  possibility 
that  F0S  < 1. 

We  see  from  (2.  6)  and  (2. 10)  that  h begins  to  increase  and  that  Fa  begins 
to  decrease.  These  behaviors  continue  toward  the  mouth  and  h can  never 
decrease  as  required  ultimately  by  physical  considerations.  We  conclude 
that  if  energy  is  conserved  along  the  estuary,  conditions  must  be  supercritical 
everywhere  and  h must  decrease  montonically. 

Another  possibility  exists  however,  namely  that  the  flow  starts  out  sub- 
critical  with  h increasing  downstream  and  then  changes  discontinuously  at  some 
section  (with  loss  of  energy  in  an  internal  hydraulic  jump)  to  a smaller  h and 
supercritical  conditions,  with  h then  decreasing  with  farther  distance  downstream. 
To  see  if  this  is  possible  we  must  make  two  investigations.  In  the  first  we 
calculate  the  momentum  balance  in  a jump  as  in  Fig.  2. 1.  We  consider  the  rate 
of  change  of  x-momentum  of  the  fluid  contained  between  Sections  A and  B at 
time  t.  At  time  dt  it  is  between  sections  A'  and  B'.  Its  rate  of  change  of 
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h.  h.-|£(h.-hb) 

W Pfg(h,-z)dz  - W [*  p,g  h,-  •^(h.-hj-z  . dz 

" 0 p,  J 

-Wj  pfg[h,-(l+ -^)zjdz 

o 1 

Integration  yields 

A.  - M.  (ha-Ksi 

g g 2Pf  'na  ' 

Using  continuity,  u,h,  = ubhb  we  may  write  this  in  either  of  two  ways: 

pHi^1  + il;  <2-14> 

F!-ltn1  + fe>  (2-,5) 

Since  hb  < h, , Ff  < 1,  >1  so  that  a jump  downstream  from  sub-critical  to 

supercritical  is  possible  from  a momentum  viewpoint. 

Let  us  now  consider  energy.  We  may  choose  a representative  sub-critical 
initial  flow  F<f  = 0. 25.  Then  if  energy  is  conserved  E'  = 2. 25  in  Eq.  (2.  8). 

Fig.  2.  2 contains  this  curve  plus  one  corresponding  to  a lower  energy  level 
E'  - 1.  50.  We  see  that  there  is  a possible  jump  to  a lower  energy  from  the 
subcritical  branch  as  shown  by  the  arrow.  We  conclude  that  a subcritical  flow 
at  the  tip  of  the  wedge  is  possible  with  h increasing  very  gradually  initially  but 
there  must  then  be  a hydraulic  jump  to  a smaller  h followed  by  a monotonic 
reduction  of  h as  the  upper  fluid  spreads  out  over  the  wedge. 


(2.12) 


(2.13) 


Fig.  2.  2 Energy  carves  for  flow  above  a salt  wedge. 
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3.  Estuary  with  a Control  Section. 


We  may  use  some  of  the  results  of  Section  2 to  discuss  an  estuary  with 
a narrows  at  the  seaward  end  as  in  Fig.  3. 1.  We  again  consider  zero  friction 
and  mixing  so  that  the  lower  fluid  will  be  inert.  In  the  region  of  the  narrows 
the  Bernouilli  equation  is  the  same  as  the  dimensional  form  of  (2.  6), 


h/A>  + — = const 


Using  continuity  uhW  = qx  = const,  we  may  write 


h&  + — — = const 
2h3  W3 

Differentiating  with  respect  to  x,  we  get 


dx‘-  h3  w3^ 


^W3  dx 


At  the  narrows  W is  a minimum  so  that  dW/dx  = 0.  If,  in  addition,  we  make  the 
reasonable  assumption  that  the  interface  slopes  upward  at  the  narrows  as  the 
issuing  fresh  water  spreads  over  the  salt  water  in  the  sea,  we  find  that  F3  = = 1 

at  the  narrows  so  that  conditions  become  critical  at  the  mouth  and  supercritical 
beyond  this.  As  we  have  seen  in  Section  2,  the  fresh-water  layer  will  continue 
to  thin  as  it  moves  seaward. 

An  interesting  relation  is  obtained  when  the  estuary  width  is  much  larger 
than  the  width  of  the  narrows.  Then  from  (3.  2) 


hcAb  + 


2Wfhf 


= DAb 
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where  the  subscript  "c"  denotes  conditions  at  the  narrows  and  D is  the  depth  of 
the  halocline  in  the  main  body  of  the  estuary.  Using  Fca  = 1 we  obtain  the  simple 
relation 

h„  =|D  (3.5) 

This  result  was  first  obtained  by  Binney  (1972). 


4.  Strongly  Mixed  Estuaries 

We  have  introduced  the  subject  of  estuaries  by  discussing  a highly  idealized 
situation  in  which  there  is  no  mixing  (and  no  friction).  Sources  of  mixing  exist  in 
all  real  estuaries,  however,  and  we  may  discuss  two  causes  in  addition  to  the 
tendency  for  breaking  waves  on  the  halocline  as  described  above. 

(1)  Influence  of  the  tides.  If  the  estuary  is  subject  to  tides,  these  will 
cause  tidal  currents  to  move  in  and  out  of  the  estuary.  As  the  current  moves 
along  the  bottom  and  sides  a turbulent  boundary  layer  will  be  created  and  there 
will  be  a tendency  to  mix  near  the  boundaries  and  cause  a reduction  of  the  density 
gradient  in  those  regions.  The  density  field  changes  momentarily  as  represented 
by  the  dashed  lines  in  Fig.  4. 1 which  portrays  a typical  pattern  with  a pronounced 
halocline.  The  new  density  distribution  will  tend  to  cause  the  flow  pattern  indicated 
by  the  arrows  in  Fig.  4. 1.  The  result  is  to  separate  the  constant  salinity  surfaces 
in  the  interior  and  to  create  a salt  flux  and  mass  flux  from  below  to  the  region  of 
the  upper  layer.  Since  we  are  dealing  with  a steady  state  over  che  long  run, 
processes  must  exist  to  restore  the  original  density  gradient  and,  in  particular, 
to  maintain  the  halocline.  The  loss  of  salt  in  the  lower  layers  is  compensated  by 
an  influx  of  salt  from  the  sea  due  to  the  estuarine  circulation  that  develops. 


z = H +AH(x  ) 


i 
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The  tides  will  also  cause  the  development  of  internal  waves.  These  may 
break  against  the  sides  of  the  estuary  again  creating  turbulence  there  and  the 
same  phenomenon  described  above  and  depicted  in  Fig.  4. 1.  Stigebrandt  (1976) 
has  conjectured  that  internal  waves  due  to  the  tides  propagate  over  the  sill  of 
the  Oslofjord  and  break  against  the  sides  and  the  head  of  the  estuary.  He  has 
constructed  an  experiment  that  demonstrates  this  effect. 

(2)  Wind  effects.  When  the  wind  blows  over  the  surface  of  the  water,  the 
resultant  agitation,  for  example  breaking  surface  waves,  will  cause  turbulence 
in  the  upper  levels.  To  the  extent  that  this  penetrates  to  lower  layers,  for 
example  to  the  halocline,  mixing  will  result.  In  addition  the  wind  will  cause 
currents  in  the  water  and  the  resulting  shear  will  be  a source  of  energy  for 
turbulence.  In  addition  the  stress  of  the  wind  may  cause  internal  waves  in  the 
deeper  layer  which  may  break  in  the  interior  or  against  the  bottom  and  sides. 

Obviously  a quantitative  representation  of  mixing  processes  will  be 
difficult  and  we  will  discuss  these  at  some  length  in  Section  5.  A great  simpli- 
fication results  if  the  mixing  is  sufficiently  strong  to  produce  a thoroughly  mixed 
estuary.  This  was  first  discussed  by  Stommel  and  Farmer  (1952,  1953).  Wc 
may  discuss  the  problem  with  reference  to  Fig.  4.  2 which  shows  conditions  in 
the  vicinity  of  the  mouth  of  an  estuary  (on  the  left)  connecting  the  estuary  to  the 
open  sea  (on  the  right).  We  have  assumed  an  estuarine  circulation  near  the  mouth 
with  fluid  of  mean  buoyancy  b,  flowing  out  with  mass  flux  qx  over  the  incoming 
fluid  with  the  buyancy  b0  of  sea  water  and  mass  flux  q0.  In  this  portion  of  the 
analysis  we  do  not  need  to  assume  uniform  buoyancies  in  the  two  layers  although 
we  will  do  this  later.  The  mean  depth  of  the  water  averaged  across  the  channel 
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is  h , the  mean  width  is  W.  Other  lengths  required  to  specify  the  geometry  of  the 

estuary  and  its  mouth  may  be  denoted  by  Lx , L There  is  a fresh-water 

flux  R into  the  estuary  and  some  kind  of  mixing.  In  the  present  discussion  the 
precise  nature  of  the  mixing  processes  is  not  important  (as  we  will  see)  but  to 
be  definite  and  simple  we  will  assume  the  mixing  is  caused  by  a rotating  propellor 
with  angular  speed  w.  Its  geometry  and  location  are  specified  by  lengths  a^aa, . . . . 

A quantity  of  importance  in  the  discussion  is  the  buoyancy  difference  b0-bj.  = 
From  dimensional  analysis  we  may  write 


= MQf.  n,  ...),  Qf  = -Att  ’ n = 

0 bo  h 2 W (hbo)* 


(4.1) 


where  a number  of  non-dimensional  ratios  of  lengths  of  the  form  I^/h , L2/h, . . . , 
al/h,  a2/h, . . . have  been  omitted.  Certainly  Ab/b0  will  vary  with  the  quantity  n; 
indeed  Ab/ho  will  decrease  as  D increases  from  the  maximum  value  of  one  when 
the  mixing  is  zero.  Ultimately,  if  the  estuary  has  finite  dimensions,  a sufficiently 
large  n , say  0. , will  cause  the  estuary  to  be  thoroughly  mixed  and  any  further 
increase  of  Q will  not  further  decrease  Ab  and  it  will  assume  its  minimum  value 
Ab..  Since  9(Ab/b0)/3Q  when  Ab  = Ab.,  we  get 


Ab. 

So 


(4.2) 


The  ratio  Tl  of  the  mean  depth  of  the  lower  layer  to  h is 


~ f a (Q t«n  ,. .. ) 


(4.3) 


When  the  estuary  is  just  thoroughly  mixed  and  Ab  = Ab. , any  further  increase 
in  the  mixing  n can  only  affect  conditions  at  the  mouth  by  the  tendency  to  further 
increase  the  turbulence  in  the  outflowing  fluid.  This  effect  is  difficult  to  analyse 
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and  we  will  neglect  it.  Subsequently  we  will  assume,  in  fact,  perfect  fluid  flow 
in  the  vicinity  of  the  mouth  so  that  the  turbulence  is  assumed  to  be  negligible 
compared  to  the  mean  motions.  Then  for  n > G.  , 71  assumes  a fixed  value  Tl. 
such  that 

\ =MQr....)  (4.4) 

In  general  we  notice  that  a combination  of  (4. 1)  and  (4.  3)  permits  us  to  write 

v.  = fs  (Q  f > 11, . • • ) (4.5) 

u0 

for  arbitrary  mixing. 

Let  us  attempt  to  find  the  form  of  the  function  in  Eq.  (4.  5)  subject  to 
several  simplifying  assumptions.  With  reference  to  Fig.  4. 2 and  Fig.  4.  3 the 
Bernoulli  equations  in  the  two  layers  are 

2 

^ + gAH  = const  (4.  6) 

2 

+ gAH  - hj.  Ab  = const  (4.  7) 

where  ht  is  the  distance  from  the  free  surface  to  the  interface,  and  where  we 
assume  irrotational  motion.  Let  us  justify  the  neglect  of  the  other  two  components 
of  the  kinetic  energy  in  (4.  6)  and  (4.  7).  The  equation  of  continuity  leads  to 


u v w 

“ »v  • 

X W h 


(4.8) 


where  X,  W 
directions. 


and  h denote  length  scales  in  the  longitudinal,  transverse 
In  each  layer  va  and  wa  will  be  negligible  compared  to  ua 


W* 


t*«l 


«1 


and  vertical 
if 


(4.9) 


We  assume  that  the  channel  has  these  properties. 


' 


Eliminating  AH  between  Eqs.  (4.  6)  and  (4.  7)  and  introducing  the  constant 
fluxes  ql  and  cjo,  we  have 


2 AbBo2  [H-hj.  -£  (x)  ]3  2 AbBx2  hf 


hx  = const 


(4. 10) 


where  B0  is  the  width  of  the  channel  at  the  level  of  the  interface,  Bj,  is  the  average 
width  in  the  upper  layer  and  £ is  the  average  height  of  the  bottom  of  the  channel 
above  z = 0.  Let  us  now  differentiate  Eq.  (4. 10)  with  respect  to  x.  We  get 

Fj-1)  =-{§  F?  + !o  d-&><H-h1-0F“  -i  Shh.F/  (4.11) 
where  F0  and  Fx  are  the  densimetric  Froude  numbers  defined  by 


*o  = — 3 » *i  = — rf 

AbB02[H-h1-C]3  toBj 


AbB^hf 


(4. 12) 


Suppose  we  have  a channel  with  a minimum  of  width  at  a certain  section 
(i.  e. , dBo/dx  = 0,  dE^/dx  = 0,  presumably  at  the  same  section).  If  the  bottom 
is  level  there  (or  if  the  mean  depth  is  also  a minimum  there),  the  left-hand  side 
of  Eq.  (4. 11)  is  zero  at  that  section.  Since  we  would  expect  a strong  slope  of  the 
interface  as  the  upper  level  moves  out  of  the  estuary  and  spreads  laterally  and 
thins  vertically  in  the  widening  channel,  we  infer  that  conditions  will  be  critical 


at  the  section,  i.  e. 


Fo  + F2  = 1 


(4. 13) 


Such  a section  is  referred  to  as  a control  jection.  Experiments  by  Stommel 
and  Farmer  (1953),  Assaf  and  Hecht  (1974)  and  Assaf,  Anati  and  Siegenthaler 
(1975)  have  shown  that  critical  conditions  do  tend  to  occur  in  a variety  of  straits 
in  laboratory  models.  If  the  strait  is  long,  friction  become  important  and  critical 
conditions  occur  at  the  ends  of  the  strait  (Assaf,  Anati  and  Siegenthaler,  1975). 


V 
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At  the  control  section  we  let  W0  and  Wx  be  the  values  of  B0  and  Bx  and  T| 
be  the  ratio  of  the  mean  thickness  of  the  lower  fluid  to  the  mean  depth  of  the  water 
(H-  C)  = h both  evaluated  at  the  control  section.  We  also  use  the  equations  for 
conservation  of  mass  and  buoyancy  (salt) 

Qo  + R = (4. 14) 

q0b0  =q1b1  (4.15) 

where  R is  the  fresh  water  discharge.  The  critical  condition  becomes 


where 


(I-ttF 


Qf  = 


R2 


boh^2 


(4. 16) 


(4.17) 


For  a given  value  of  Qf , Eq.  (4. 16)  represents  a curve  like  that  shown  in 
Fig.  4.4.  It  is  of  the  form  of  Eq.  (4.  5)  yielding  £b/b0  as  a function  of  T]  instead 
of  the  mixing  parameter  n.  Physically,  if  we  increase  0 from  a small  value, 
we  expect  Ab/b0  to  decrease  monotonically  from  the  value  1 corresponding  to 
zero  mixing.  T1  will  increase  or  decrease  with  increase  of  Q , and  we  move 
on  the  curve  of  Fig.  4. 4 toward  the  minimum  point.  At  a certain  value  n.  the 
estuary  will  be  thoroughly  mixed,  £b  will  no  longer  decrease  and  9(Ab/b0)/9Q  = 0. 
We  have 


3(Ab/bn)  8(Ab/bo)  97] 

an  ~ ail  an 


(4.18) 


Obviously  if  9(£b/b0)/971  = 0,  corresponding  to  the  minimum  point  of  the  curve 
in  Fig.  4.4,  9(Ab/b0)/9Q  = 0 and  Ab/b0  will  have  the  value  appropriate  to  the 
thoroughly  mixed  estuary.  The  condition 
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_0 

an  ~u 


(4. 19) 


together  with  Eq.  (4. 16)  determines  the  problem  of  a thoroughly  mixed,  or  in 
the  terminology  of  Stommel  and  Farmer,  an  "overmixed"  estuary.  Notice  that 
the  conservation  equations  (4.14),  (4.15)  yield 


UU  — 

<!‘br  = R 


(4.  20) 


so  that  a minimum  of  Ab/b0  corresponds  to  a maximum  of  the  discharge  qx . 

The  variation  of  n2  with  n depends  on  the  particular  geometry  of  the 
estuary,  and  although  this  variation  may  have  some  quantitative  importance  we 
will  consider  instead  the  idealized  channel  with  vertical  sides  so  that  n*3  =1. 


Eq.  (4. 19)  then  yields 


i- = X 
v b0/  (i-n) 


In  terms  of  the  densimetric  Froude  numbers  we  may  write 

*o  = n , f?  = i - n , q r = 


(4.21) 


(4.  22) 


One  special  case  involves  the  condition  in  which  the  densities  of  the  two  fluids  are 
nearly  equal.  Then  from  (4.  21),  Tl  ~ §,  R~0.  Experiments  by  Stommel  and 
Farmer  (1953)  and  Assaf , Anali  and  Siegenthaler  (1977)  confirm  the  prediction 
that  the  two  fluids  have  equal  thicknesses  in  thoroughly  mixed  estuaries  with  small 
fresh  water  influx  and  I have  derived  this  result  theoretically  (Long,  1976).  A 
second  special  case  is  one  in  which  T)  -•  0 so  that  no  salt  water  enters  and  fresh 
water  fills  the  estuary.  In  this  case  the  solution  yields  Ab/b0  =1,  Fq  = 0 as 
also  required  by  physical  considerations.  In  the  general  case  (4.  22)  shows  that 
Tl  < i so  that  the  interface  is  always  at  or  below  mid-depth. 
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Of  interest,  particularly  in  later  discussions,  is  the  height  AH©  of  the 
water  surface  in  the  estuary  above  the  level  of  the  ocean  surface.  If  we 
evaluate  the  constants  in  Eqs.  (4.  6)  and  (4.  7)  in  the  estuary  and  in  the  ocean 
respectively,  we  get 


f + gAH  = gAHo 
2 

+ gAH  -Ahi  = 0 


(4.  23) 


(4.  24) 


where  we  now  assume  that  the  widths  in  the  ocean  and  estuary  are  very  large 
compared  to  the  width  of  the  channel.  The  constant  in  Eq.  (4. 10)  is  -gAH0/&> 
and  this  equation  may  then  be  written 


where 


Using  (4.  22)  we  get 


Sfa-  ^ - <1-1  > + B = 0 


B = — 

hAf/pf 


3 = il-y'n  j (4.27) 

In  the  overmixed  state  6 ranges  between  i and  f. 

We  may  contrast  this  strongly  mixed  case  with  the  case  considered  in 
Section  3 in  which  the  mixing  is  zero.  Then  F§  = 0,  F*  = 1 and  Eq.  (4. 25) 
yields  8 = •§  (1-T| ).  For  future  reference  we  list  the  two  behaviors: 


(4.25) 


(4.  26) 


Fo  =°*  Pf  =1,  3 (1-H),  Qfa  = (l-Tl)3 


= 1 . = 1 , weak  mixing 


(4.28) 


( 


F<?  = n . f l = a-n ).  e = la-  §■  n >.  o?  = [p^3 


qi  _ (lili?  Ab  _ (i-2jij 

R (1-2T1 ) ’ b0  (1  — T)  )*  • Str°ng  miXing 


(4. 29) 


5.  Turbulence-  and  Mixing  Processes  in  Stratified  Fluids. 

We  have  decided  that  turbulent  mixing  of  waters  of  different  density  (salinity) 
is  a fundamental  process  in  the  dynamics  of  estuaries  and  we  now  consider  some 
aspects  of  turbulence  in  stably  stratified  fluids. 

Turbulence  in  a fluid  needs  an  energy  source  or  it  will  die  out.  This  is 
especially  true  in  a stably  stratified  fluid.  Thus,  letting  primes  denote  departures 
from  mean  conditions  and  using  the  equations  of  motion,  we  may  form  the  energy 
equation  for  the  turbulent  kinetic  energy 

T'  = .u'a+v^a+w'3  (5.1) 

If  conditions  are  statistically  steady  and  horizontally  homogeneous,  we  get 

Jf(T')  ="  W(T'+p'/pf)]  - Wu,  ->Pb'  - c (5.2) 

where  e is  the  dissipation  function,  p is  pressure  and  subscript  "z"  denotes  the 
vertical  derivative.  When  the  fluid  is  stably  stratified,  i.  e. , if  the  mean  density  (or 
buoyancy)  decreases  with  height,  the  turbulence  will  be  accompanied  by  a 
buoyancy  flux  q = -w'b'  < 0 so  that  this,  and  the  viscous  dissipation  e,  act  as 
sinks  of  kinetic  energy.  Energy  must  come  from  the  first  and/or  the  second  term 
on  the  rhs  of  Eq.  (5.  2).  The  first  is  the  energy-flux-divergence  term  and  is  a 
basic  source  term  in  certain  laboratory  experiments  to  be  discussed  in  this  section. 
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The  second  is  always  a source  term  and  arises  from  a conversion  of  energy 
of  the  mean  shear  of  the  current  to  turbulent  energy.  It  is  believed  to  be  of 
fundamental  importance  in  geophysical  phenomena. 

In  a stably  stratified  estuary  in  which  salinity  increases  with  depth , q 
will  tend  to  be  negative  because  rising  parcels  tend  to  be  more  salty  and  heavy 
and  falling  parcels  tend  to  be  less  salty  and  lighter.  One  must  be  careful, 
however,  because  of  the  possibility  of  wave  motions  contributing  to  w'  and  b'. 
Thus,  if  the  fluid  is  stably  stratified,  it  is  capable  of  internal  gravity-wave 
motion  in  which  the  basically  level  density  surfaces  move  up  and  down  in  waves. 
Obviously  if  the  waves  do  not  break,  there  will  be  no  rupture  of  these  surfaces 
and  therefore,  neglecting  the  very  small  molecular  conduction,  no  flux  of 
buoyancy  (or  salt)  despite  sizable  values  of  w'  and  b'.  The  correlation  co- 
efficient will  be  zero.  If  the  v'aves  break,  there  will  be  intermittent  turbulence 
superimposed  on  the  wave  motion  and  q will  be  negative  although  the  correlation 
coefficient  may  be  much  less  than  one.  Negative  q means  that  the  kinetic 
energy  tends  to  decrease  because  it  requires  work  to  lift  heavy  parcels  up 
and  bring  light  parcels  down.  There  is  a tendency  in  doing  this  to  increase 
potential  energy  at  the  expense  of  kinetic  energy. 

To  demonstrate  this,  it  is  useful  to  define  available  potential  energy  per 
unit  mass  (Long,  1970)  by  considering  it  to  be  the  kinetic  energy  per  unit  mass 
attained  by  a parcel  of  buoyancy  b = b'+b(z)  as  it  falls  from  the  height  z to  the 
height  z0  at  which  its  buoyancy  b is  equal  to  the  mean  buoyancy  b(z0)  at  that 
level.  We  have  b'  =b(z0)-b(z)  =-bzf  approximately,  where?  = z-z0  and  we 
have  assumed  that  ? is  small  compared  with  the  length  scale  of  the  vertical 
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variation  of  mean  buoyancy1 * . Then  neglecting  disturbance  pressure,  we  may 
write 


dw . d3e 
dt  ~ "dt* 


= -b'  =bI? 


(5.3) 


because  z0  is  a Lagrangian  quantity,  so  that  dzo/dt  = 0. 


Integrating,  we  get 


w3/2  - b,e3/2  = const  (5.4) 

Thus  available  potential  energy  may  be  defined  as 

V'  = -htf/Z  or  V’  = b’5/2  (5.  5) 


We  may  also  identify  q with  potential  energy  changes.  The  potential  energy 
of  a particle  of  volume  V0  and  density  p is  PgV0z.  Let  us  now  define  incremental 
potential  energy  as  pgV0z-pfgV0z  so  that  this  potential  energy  is  zero  when  the 
particle  has  the  density  p{.  If  we  let  V represent  the  incremental  potential 
energy  per  unit  mass,  then  to  within  the  Boussinesq  approximation,  V = bz. 

Since  b is  nearly  conservative,  putting  b = b'+  b and  assuming  no  mean  vertical 
velocity,  we  have 

dV/dt  = wV  = -q  (5.  6) 

is  the  average  rate  of  increase  of  incremental  potential  energy  per  unit  mass. 

We  may  identify  this  with  available  potential  energy  by  differentiating  the  second 
equation  in  (5.  5)  and  again  assuming  db/dt  = 0.  We  get 

dV’/dt  =ib'd?/dt  - i^db/dt  = \ w'b'-i^bjW'  = w'b'  (5.7) 

so  that 

dV'dt  = w'b'  =-q  (5.8) 

1 This  may  not  always,  or  even  usually,  be  the  case  but  our  development  here 

is  only  suggestive  of  the  definition  in  Eq.  (5.  5). 
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Comparing  (5.  6)  and  (5.  8),  we  see  that  the  average  rate  of  increase  of  incremental 
potential  energy  and  the  average  rate  of  increase  of  available  potential  energy  are 
the  same. 

In  application,  we  can  conceive  of  an  energy-containing  eddy  in  the  form  of 
a whirl  with  horizontal  axis  of  rotation,  with  velocity  cr„  and  diameter  l.  It  will 
lift  parcels  from  their  level  of  origin  a distance  5 ~ X so  that  V'  ~ <7b  X , where 
crb  is  the  rms  buoyancy  fluctuation.  If  the  turbulence  is  not  decaying,  the  kinetic 
energy  must  be  of  this  order  or  larger  so  that  ob  X/o*  l. 

A typical  phenomenon  in  a turbulent,  stably  stratified  fluid  is  the  appearance 
of  thin  layers  across  which  density  changes  abruptly.  We  have  seen  that  these 
exist  in  nature,  for  example  the  interface  between  the  river  water  and  the  salt 
water  in  a salt  wedge  at  the  mouth  of  some  rivers.  The  surprising  feature, 
however,  is  that  turbulence  actually  sharpens  the  interface,  or  at  least  serves 
to  maintain  the  interface  instead  of  diffusing  it  as  one  might  expect  from  the 
general  diffusive  nature  of  turbulence.  An  important  example  is  an  experiment 
originally  by  Rouse  and  Dodu  (1955)  portrayed  in  Fig.  5. 1.  The  vessel  is  filled 
initially  with  water  with  a stable  linear  density  profile.  The  grid  is  then  activated 
and  turbulence  is  created  in  the  nearby  fluid.  The  result  is  a growing  turbulent 
upper  layer  with  a nearly  constant  mean  density  separated  from  the  quiescent 
fluid  below  by  a thin  interfacial  layer  of  large  density  gradient.  We  sometimes 
idealize  this  by  considering  it  a density  discontinuity.  Actually  in  this  experiment 
the  layer  has  a thickness  of  1-2  cm.  In  general,  the  thickness  is  proportional  to 
the  depth  of  the  mixed  layer  with  a constant  of  proportionality  of  1/5-1/6.  (Moore 
and  Long,  1971,  Wolanski,  1972,  Long,  1973,  Crapper  and  Linden,  1974,  Wolanski 
and  Brush,  1975,  Assaf,  Anati  and  Siegenthaler,  1977)  The  interface  moves  away 
from  the  grid  at  a speed  called  the  entrainment  velocity  u, . 
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A detailed  understanding  of  the  turbulence  in  this  experiment  has  not  yet 
been  achieved  but  Linden  (1973)  has  performed  allied  experiments  and  has 
suggested  that  the  large  eddies  in  the  upper  mixed  layer  deflect  the  interface 
downward,  storing  potential  energy.  When  this  is  released  by  upward  motion, 
a portion  of  the  heavier  fluid  is  ejected  into  the  mixed  layer  and  then  carried 
away  by  the  turbulent  eddies,  or  as  in  Keulegan's  experiment  (1949)  by  the  mean 
current,  leaving  the  interface  sharp  again. 

The  experiment  of  Rouse  and  Dodu  is  typical  of  those  without  shear  and  the 
only  source  of  energy  is  the  energy-flux-divergence  term  of  Eq.  (5. 2).  We  will 
see  that  similar  experiments  have  been  constructed  with  shear  and  in  these  the 
shear  term  of  Eq.  (5.  2)  is  an  important  energy  source. 

Cromwell  (1960)  studied  an  experim  3nt  similar  to  that  of  Rouse  and  Dodu 
to  simulate  the  pycnocline,  but  the  first  reliable  data  were  obtained  by  Turner 
(1968)  in  an  experiment  in  which  the  lower  fluid  was  agitated  and  fluid  was  with- 
drawn from  the  stirred  layer  at  a rate  adjusted  to  keep  the  interface  at  the  same 
distance  from  the  grid.  The  entrainment  velocity  is  then  defined  by  Au,  = Q 
where  Q is  the  volume  withdrawn  per  unit  time  and  A is  the  cross-sectional 
area  of  the  tank. 

There  have  been  a number  of  recent  experiments  similar  to  those  of 
Rouse  and  Dodu  and  of  Turner,  for  example  by  Brush  (1970).  Equipment 
identical  to  that  of  Turner  was  constructed  by  Wolanski  (1972)  (see  also  Wo’  - nski 
and  Brush,  1975),  and  the  one-and  two-grid  experiments  were  run  with  strati- 
fication caused  by  heat,  salt,  sugar,  suspensions  of  sediments  and  minute  silica 
spheres.  Additional  experiments  have  been  run  in  Turner's  apparatus  by  Linden 
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(1973),  Crapper  (1973)  and  Crapper  and  Linden  (1974)  using  heat  and  salt.  Baines 
(1975)  has  studied  a similar  experiment  with  entrainment  caused  by  a jet  impinging 
on  an  interface. 

An  important  finding  in  the  experiments  by  Turner  (1968)  may  be  expressed 

as 

U8/u>  =C[a/7(Ab)]n  (5.9) 

where  u>  is  the  frequency  of  the  oscillating  grid  and  C is  independent  of  u>  and  £b. 

A number  of  lengths  are  kept  constant  in  the  experiment  so  that  the  dimensional 
quantity  C may  vary  with  these.  Turner  found  that  for  larger  values  of  Ab/w*5, 
the  exponent  n = 1 when  stratification  is  caused  by  temperature  differences  and 
n = 3/2  when  caused  by  differences  in  salt  content.  Later  investigations  (Wolanski, 
1972)  have  confirmed  these  results  and,  very  recently,  Crapper  and  Linden  (1974)  have 
shown  rather  convincingly  that  the  difference  in  the  values  of  n is  due  to  the  influence 
of  the  relatively  large  molecular  conductivity  kh  in  the  heating  experiments  (the 
coef'icient  ks  is  much  smaller  for  salt).  It  appears  that  whenever  the  n = 1 law 
describes  the  entrainment,  the  thin  layer  of  strong  density  variation  has  an  inner 
layer  or  core  in  which  molecular  diffusion  is  important.  Indeed,  earlier  unpublished 
experiments  by  Claes  Rooth  (Turner,  1973)  support  this  interpretation.  Rooth 
found  a 3/2  dependence  in  heating  experiments  when  larger  turbulent  velocities 
are  generated.  Thus  it  appears  well  established  that  the  3/2  dependence  is  appro- 
priate for  larger  Peclet  numbers,  Pe  = ~u  Vkh  or  <7ttt/k,  where  and  i arc  the 
velocity  and  length  units  of  the  turbulence.  Crapper  and  Linden  suggest  a thresh- 
hold  value  of  Pe  ~ 200  when  = <j'a  and  H = V are  characteristic  of  the  turbulence 
near  the  interface.  The  dependence  on  Reynolds  number.  Re,  has  net  been 
established  because  of  the  rather  small  ranges  of  Re  in  the  experiments,  but  both 
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Wolanski  (1972),  who  varied  Re  by  a factor  of  threes  and  Crapper  and  Linden  (1974) 
report  very  weak  dependence,  if  any1 . We  may,  therefore,  write  for  large  Pe  and 
Re,  and  strong  stability, 

3 

u./«  = Cxw3  / (£b) 2 

where  Cx  is  a function  of  a,  D,  and  the  lengths  aj, , • • » characteristic  of  the 

grid.  It  is  convenient  to  introduce  a dimensionless  quantity  K„  by  the  definition 

C4  = a4D-2K„  (a/D,a/ai,a/a2. . ) (5.10) 

and,  therefore, 

u,/u*  =KnRi+'?,  Ri*=DAb/u*  (5.11) 

where  u*  = u>a.  It  is  likely  that  K„  is  independent  of  a/D  when  this  ratio  is  small. 

Ri  is  called  the  overall  Richardson  number. 

The  Turner  or  Rouse  and  Dodu  experiment  can  be  run  with  an  initial  two-fluid 
system  in  which  both  layers  have  a nearly  uniform  but  different  density.  The  inter- 
face again  moves  away  from  the  grid  and  the  system  is  basically  unsteady  because 
the  density  of  the  upper  fluid  is  increasing  and  because  the  interface  is  moving. 

The  motion  of  the  interface  can  be  eliminated  by  introducing  fluid  of  the  density  of 
the  lower  layer  into  the  lower  layer  at  a rate  q0  = Aue . The  upper  surface  of  the 
system  can  be  held  stationary  by  removing  an  equal  volume  of  the  upper  brackish 
water.  There  will  still  be  an  unsteadiness,  however,  because  the  density  of  the 
upper  fluid  will  still  increase.  This  can  be  eliminated  by  adding  fresh  water  to  the 
upper  layer  with  a flux  R and  removing  brackish  water  at  the  rate  qa . The  following 


1The  independence  of  molecular  quantities  is  common  in  turbulence  (Tennekes 

and  Lumley,  1972).  In  Turner's  original  paper  (1968),  he  expressed  the  belief  that 
the  n = 1 law  was  the  fundamental  one  and  that  in  some  manner  the  very  low  diffusivity 

of  salt  caused  the  n = 3/2  law.  He  has  changed  his  mind  on  the  basis  of  the  evidence 

we  present  here  (Turner,  1973). 
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relations  must  be  satisfied  to  insure  a steady-state 


q0  =Au. 

(5. 12) 

q0+R  = qx 

(5. 13) 

Qo^o  =Qil>i 

(5. 14) 

An  analogous  experiment  was  set  up  by  Welander  (1974).  It  consisted  of 
a large  vessel  (Fig.  5.  2)  with  a small  semi-enclosed  portion.  The  vessel  was 
filled  with  salt  water  initially.  Fresh  water  was  then  introduced  into  the  semi- 
enclosed  region  and  the  upper  portions  of  this  region  was  stirred  with  a propcllor. 
A steady  state  developed  with  an  interface  in  the  semi-enclosed  region  at  a depth  I) 
and  with  all  conditions  in  Eqs.  (5. 12)-(5.14)  automatically  satisfied.  The  system 
resembles  an  estuary -ocean  system;  indeed,  this  was  the  basic  objective  of  the 
experiment.  Welander  only  measured  the  depth  D as  a function  of  R and  we  will 
discuss  these  observations  in  Section  6. 

When  the  interface  is  stationary,  as  in  an  estuary  or  in  Turner's  or 
Welander's  experiment,  there  is  a flow  of  mass  through  the  interface  with  a 
mean  vertical  velocity  u, . We  may  then  consider  the  buoyancy  budget  in  an 
element  of  the  thickness  of  the  interfacial  layer,  arbitrary  length  and  width  of  the 
estuary  or  vessel  (Fig.  5.3).  Integrating 

Jjp  + v.  v b = 0 

over  the  element,  assuming  steady  state  and  horizontal  homogeneity,  we  get 

-(wb)o  + (wb)i  - (w'b')o  + (wT>')x  = 0 

Since  the  element  is  very  thin,  we  may  neglect  the  variation  of  the  mean  vertical 
velocity  from  the  lower  to  the  upper  surface  and  set  w = u,  . Also  the  lower  layer 
is  non-turbulent  so  that  (w'b')0  = 0.  We  get 
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q = -u.  tt) 


(5.15) 


where  q is  the  buoyancy  flux  just  above  the  interface. 

Several  experiments  have  been  constructed  to  inti  oduce  shearing  currents 
in  turbulent  density- stratified  systems  in  an  effort  to  simulate  atmospheric  and 
oceanic  phenomena.  The  first  of  these  of  direbt  relevance  to  our  discussion  was 
that  of  Kato  and  Phillips  (1969).  The  apparatus  was  a large  circular  annular 
channel  filled  with  salt  water  with  an  initial  linear  density  gradient.  A constant 


stress  t = u*  was  applied  by  rotating  a flat  screen  at  the  surface  (Fig.  5.  4).  For 


larger  values  of  Ri  , they  found 


U./U*  = K„Ri 


(5.16) 


where  Ri  is  of  the  same  form  as  in  Eq.  (5. 11)  and  is  the  buoyancy  jump 
from  the  upper  mixed  layer  to  the  quiescent  region  below1 . 

An  experiment  by  Moore  and  Long  (1971)  was  constructed  to  permit  a 
steady  state.  In  a large  channel  shaped  like  a race  track,  fluid  was  injected 
from  nearly  horizontal  jets  at  bottom  (salt  water)  and  top  (fresh  water)  in  opposite 
directions  to  obtain  a shearing  current  (Fig.  5.  5).  Mean  zero  vertical  velocities 
were  achieved  by  withdrawing  equal  volumes  of  fluid  through  numerous  holes  at 
bottom  and  top.  At  larger  values  of  the  density  difference,  two  homogeneous 
layers  existed  at  top  and  bottom  with  an  interface  in  the  middle.  The  salt  water 
in  the  jets  comes  from  a reservoir  and  the  withdrawn  fluid  at  the  bottom  is  pumped 
back  into  the  reservoir  which  is  kept  at  a constant  level.  The  jets  at  the  top  are 


1 More  recent  experiments  have  been  run  in  the  Kato  and  Phillips  apparatus 
using  a lower  fluid  of  uniform  density  (Kantha,  1975).  The  results  do  not  yield  a 
simple  power  law  but  rather  a faster  and  faster  decrease  of  u#  with  Ri*.  It  is 
possible,  however,  that  molecular  viscosity  becomes  more  and  more  important 
at  the  higher  values  of  Ri*. 


. 
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of  tap  water  and  the  slightly  salty  withdrawn  fluid  at  the  top  is  pumped  to  waste. 
Since  the  fluid  returned  from  the  bottom  to  the  reservoir  is  somewhat  less  salty 
than  that  in  the  lower  jets,  salt  must  be  added  continually  to  keep  the  reservoir 
at  a fixed  density.  The  added  salt  is  transported  vertically  by  the  turbulence. 
The  salt  flux  is  known,  of  course,  and  this  can  be  used  directly  to  compute  the 
buoyancy  flux.  The  experiment  yielded 


q = -K3(Au)3/D  (5.17) 

over  a range  of  D£b/(Au)a  of  1 - 60,  where  A u is  the  velocity  difference  between 
mean  velocities  measured  near  the  top  and  bottom  and  K3  is  a positive  constant 
of  proportionality.  If  we  define  the  entrainment  velocity  by  -14  Bo  = q,  Eq.  (5. 17) 
yields  the  same  result  as  in  Kato  and  Phillips  (Eq.  5. 16)  if,  as  seems  very  likely 
from  a discussion  of  the  Kato  and  Phillips  experiment  by  Long  (1975  a),  Au/u*  is 
independent  of  the  Richardson  number,  where  u*  is  the  constant  momentum  flux 
in  the  tank. 

Finally,  in  a recent  experiment  by  Wu  (1973),  the  source  of  energy  and 
shear  was  a current  of  air  blowing  over  a vessel  containing  a two- fluid  system 
(Fig.  5.  6).  Wu  also  obtained  Eq.  (5. 16)  although  his  coefficient  of  proportionality 
was  much  smaller.  He  conjectured  that  this  was  because  of  the  very  different  shear 
produced  in  a closed  container,  but  Bo  Pedersen  (private  communication)  has 
suggested  that  the  flow  at  the  interface  may  have  been  laminar  over  much  of 
the  length. 


ajc 

The  different  dependence  on  Ri  for  the  two  experiments  has  been  the  source 
of  perplexity  (Turner,  1973,  Linden,  1973)  because  the  mixing  processes  appear 
to  be  very  similar.  Indeed,  Linden  has  stated  that  tne  Kato  and  Phillips  data  are 
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also  consistent  with  a -3/2  behavior,  although  support  for  the  statement  seems 
lacking.  We  attempt  below  to  construct  a unified  understanding  of  the  two  results. 

Turner  (1973)  has  made  the  valuable  suggestion  that  the  erosion  of  the  inter- 
face should  depend  on  the  properties  of  the  turbulence  near  the  interface,  in 
particular  on  the  rms  velocity  oj,  and  the  integral  length  scale  V near  the 
interface.  Thus,  he  proposed  the  form 

u ./<x;  = FfRii),  m1  = (5.18) 

where  possible  dependence  on  other  quantities  is  suppressed  and  it  is  assumed 
that  Pe  and  Re  are  large.  In  an  attempt  to  determine  the  dependence  on  Rij  in 
his  density-interface  experiments,  in  which  o,,'  and  V were  not  measured,  Turner 
used  experimental  data  by  Thompson  and  Turner  (1975)  in  Turner's  apparatus 
with  a homogeneous  fluid  and  one  grid.  They  measured  cr„  and  i at  many  levels 
and  found  that  au'  was  proportional  to  w and  that  j l increased  linearly  with  distance 
from  the  grid  but  was  independent  of  w.  Although  Thompson  and  Turner's  experi- 
ment had  no  density  variation,  Turner  (1973),  Thorpe  (1973)  and  Crapper  and 
Linden  (1974)  have  assumed  that  the  results  are  directly  applicable  to  the  mixing 
experiments.  Thus,  at  z = D they  use 

crj/wa  = C3(a/D,  a/ai  ,a/aa. . . ) (5.19) 

l'/D  = C3  (a/ax , a/as . . . ) (5. 20) 

so  that  Eq.  (5. 11)  may  be  written 

ue  /or;  = K*  Ri~ K K4  = K4(a/D,  a/ai , a/a2. . . ) (5.  21) 

Neglecting  viscosity,  the  proportionality  of  o;  and  u>  follows  from  dimensional 
analysis  but  only  when  the  fluid  is  homogeneous  because  the  presence  of  an  inter- 
face introduces  a new  quantity  involving  time,  namely  Gb. 


We  may  also  obtain  a dependence  on  R ij,  for  the  shearing  experiments.  With 
shear  and  horizontal  homogeneity  the  averaged  x-equation  of  motion  yields 

9t/9z  = da/dt  (5.  22) 

where  u is  the  mean  horizontal  velocity  at  depth  z.  In  the  steady-state  experiments 
of  Moore  and  Long  (1971),  8t/9z  = 0 and  therefore  t is  constant  with  height.  Since 

t = -u'w',  and  since  the  correlation  coefficient  is  very  likely  to  be  of  order  one  in 

l 

the  mixed  layers,  it  follows  that  u+  ~ t is  proportional  to  au'.  The  interface 
introduces  the  length  D and  it  seems  reasonable  that  the  eddies  fill  the  whole 
depth  since  the  density  gradient  in  the  mixed  layer  is  weak1 . Thus  we  use  V ~ D 
and  obtain  from  Eq.  (5. 16) 

u„  /oj  = KeRij-1  (5.23) 

for  the  Moore  and  Long  experiment.  In  the  experiment  of  Kato  and  Phillips,  we 
may  use  Eq.  (5.  22)  to  obtain  the  increment  of  t over  the  depth  D.  It  is 

At/ t ~ UD/T  jU*  (5.24) 

where  U is  the  speed  of  the  screen  and  where  Tj  is  the  time  period  for  a change 
of  depth  of  order  D so  that  Tt  ~D/ue . Therefore,  ignoring  sidewall  effects, 

At/t  = (U/u*)  (ue  /u*)  ~ Ri*-1  (5.  25) 

)|( 

since,  as  shown  by  Long  (1975a)  U/u*  is  independent  of  Ri  . This  reveals  that 
the  stress  varies  very  little  over  the  depth,  that  u*  ~ and  that  Eq.  (5.  23) 
again  holds. 

1 If  a fluid  is  homogeneous,  the  large,  energy-containing  eddies  tend  to  be  as 
large  as  the  dimensions  of  the  region  (Tennekes  and  Lumley,  1972).  Here  the 
available  potential  energy  may  be  as  large  as  the  kinetic  energy  but  this  should 
not  change  the  order  of  magnitude  of  thj  eddy  size.  Indeed  visual  observations 
indicated  that  the  eddies  filled  the  whole  mixed  layer  in  the  Moore  and  Long  experiment. 
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Thus,  two  different  entrainment  velocities,  Eqs.  (5.21)  and  (5. 23),  are 
indicated  in  the  two  cases  even  when  the  characteristics  of  the  eroding  eddies 

are  the  same  and  this  is  more  perplexing  than  the  difference  in  the  exponent 

* * 
of  Ri  . However,  a different  dependence  of  ctu'/u*  or  a^/ua  on  Ri  in  experi- 
ments without  shear  may  be  obtained  by  a plausible  argument  which  casts  doubt 
on  the  applicability  of  Thompson  and  Turner's  experiment,  in  particular  Eq. 

(5. 19),  to  an  experiment  with  a density  interface.  When  there  is  shear,  we 
have  seen  that  experiment  indicates 

q ~ ojf/D  ~ u* /D  (5.26) 

Let  us  now  evaluate  q in  the  mixed  layer  near  the  interface.  We  get  q ~ ajaj 
where  ar b is  the  rms  buoyancy  fluctuation  near  the  interface,  and  we  make  the 
plausible  assumption  that  the  correlation  coefficient  is  of  order  one.  Thus 

oB'a/ay3~l  (5.27) 

Assuming  again  that  the  eddies  fill  the  whole  layer,  kinetic  energy  and  available 
potential  energy,  ob'D,  are  of  the  same  order  in  the  mixed  layer.  Although  the 
layer  has  very  little  density  variation,  this  result  may  be  obtained  by  considering 
first  an  experiment  with  very  strong  turbulence  imposed  externally.  Then 
<ra'a/olJl'  will  be  very  large.  As  we  decrease  the  turbulence  in  successive  experi- 
ments, this  ratio  will  decrease.  If  turbulence  continues  to  exist,  the  ratio  has  a 
lower  limit  because  <ju’2 /o'hV  < l would  imply  that  T'  < V'  and  a consequent 
cessation  of  the  turbulence.  It  seems  unlikely  that  the  ratio  will  get  large  again 
as  stability  increases  so  that  cru'2/cr^  should  approach  a constant.  The  argument 
is  equally  valid  with  or  without  shear.  Using  Eq.  (5. 27)  when  shear  is  absent 
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to  eliminate  in  the  relationship  q ~u,  ~ ayb  , the  -3/2P.tw  leads  to 

u,  /u * ~ <xj3  /Du„  tb  ~ u3 /(Ddb) 2 (5. 28) 


or 


CTu'/u+  ~Ri*“  B’t  u*  = 


coa 


(5. 29) 


instead  of  Eq.  (5. 19).  The  decrease  of  rms  velocity  with  increase  of  Richardson 
number  when  density  variations  are  present  may  be  caused  by  the  weak  density 
gradient  in  the  mixed  layer.  In  a layer  as  a whole,  the  slight  density  variation 
still  has  dynamic  importance  as  indicated  by  the  proportionality  of  kinetic 
energy  and  available  potential  energy  and  by  the  fact  that  q,  varying  linearly  in 
the  layer,  has  a relatively  large  value  near  the  interface.  Such  arguments  have 
been  advanced  earlier  by  the  author  (Long,  1972,  1973). 

The  energy  argument  may  be  amplified.  Rouse  and  Dodu  (1955)  and  others 
(Kato  and  Phillips,  1969,  Turner,  1973,  Wu,  1973)  have  suggested  that  the 

.*-i 

Ri  law  implies  that  the  change  of  potential  energy  is  proportional  to  the 
energy  supply  by  the  external  source.  Thus,  as  we  have  seen,  the  average 
rate  of  increase  of  potential  energy  per  unit  mass  is  q,  so  that  the  rate  of  in- 
crease of  potential  energy  for  the  system  is  proportional  to  qD.  In  the  Kato  and 
Phillips  experiment,  for  example,  the  rate  of  working  of  the  external  force 
is  tU.  If  these  are  proportional, 

q ~ tU/D  ~U^/D  (5.  30) 

as  in  Eq.  (5. 26).  The  same  conclusion  cannot  be  reached  for  cases  without 

*__.3 

shear  and  on  this  basis  it  may  be  argued  that  the  Ri  2 law  does  not  conform 
to  any  simple  energy  argument.  We  may  show,  however,  that  the  last  conclusion 
is  not  correctly  drawn.  In  the  shearing  experiments,  the  velocity  difference  is 
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proportional  to  u+  = -(w'u1)^  (Long,  1975a)  and  the  two  energy-source  terms  in 
Eq.  (5. 2)  as  well  as  the  dissipation  e,  are  of  order  of  /D  or  cr?/l'  near  the 
interface.  If  q ~ u*  Ab  is  also  of  this  order,  we  obtain  u,  /<r'u  ~ Ri^1  as  in 
Eq.  (5.  23).  When  shear  is  absent,  the  single  source  term  is  the  first  term 
on  the  rhs  of  Eq.  (5.2)  and  is  also  of  order  a? /V.  The  Ri^1  law  again  implies 
equality  of  all  sink  and  source  terms.  The  correct  interpretation  of  experi- 
mental results  thus  seems  to  be  that  the  turbulence  has  a character  that  causes 
potential  energy  to  increase  at  a rate  proportional  to  the  rate  at  which  kinetic 
energy  is  supplied  to  the  region  of  the  interface  and  not  necessarily  proportional 
to  the  rate  of  generation  of  kinetic  energy  at  the  external  source. 

An  additional  piece  of  information  may  be  added  in  relation  to  experiments 
without  shear.  If  we  assume  that  the  small  buoyancy  difference  Ab  across  the 
mixed  layer  is  of  order  of  the  rms  buoyancy  fluctuation  (implying  an  eddy  length 
scale  of  the  order  of  the  depth  of  the  layer),  Eq.  (5. 27)  and  Eq.  (5.  29)  lead  to 

Ab/Ab~Ri*~T  (5.31) 

This  quantity  was  measured  by  Wolanski  (1972)  and  Wolanski  and  Brush  (1975) 
for  the  salt  experiments  (Fig.  5.  7).  There  is  good  agreement  with  Eq.  (5.  31) 
especially  at  higher  values  of  Ri  . 

When  there  is  shear,  oi  in  Eq.  (5.  27)  is  of  order  u+  and  <r{,  is  of  order 
lb  so  that  AbD/u*  ~ 1 and  Ab/ Ab  is  proportional  to  Ri*-1 . Notice  also  that  in 
the  mixed  layer  the  velocity  shear  is  u,  ~u*/D,  as  shown  by  Long  (1975a).  Thus 
the  gradient  Richardson  number  Ri  = bt/(ut)3  is  of  order  one  in  the  mixed  layer. 
Observations  in  the  lower  mixed  layer  in  the  atmosphere  (Businger  ct_al.  1971) 
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indicate  that  this  has  a maximum  of  Rie  = 0.  20  and  we  may  speculate  that  this 
is  the  magnitude  of  Ri  in  the  mixed  layer  in  the  experiment.  As  we  have  already 
suggested,  Ri  ~ 1 means  that  the  slight  density  gradient  in  the  mixed  layer  has 
dynamical  significance. 

We  have  suggested  that 

u./<ru  =Knof/D£b  (5.32) 

is  a universal  law  for  entrainment  with  or  without  shear,  where  ott  is  the  rms 
velocity  at  or  near  the  interface.  Let  us  apply  this  to  an  experiment  without 
shear,  as  in  Turner's  experiment  except  that  the  lower  quiescent  fluid  has  a 
linear  density  gradient.  If  we  also  adopt  the  findings  of  Bouvard  and  Dumas 
(1967)  and  Thompson  and  Turner  (1975)  that  cru  is  proportional  to  D , then 
since  db  - if  D/2,  where  N2  is  the  buoyancy  gradient  in  the  lower  fluid,  we 
find  that  D <*  t2  ^ e.  This  behavior  was  proposed  by  Linden  (1975)  and  his 
experiments  with  this  type  of  fluid  system  provided  close  verification  of  this 
time  dependence.  These  experiments,  therefore,  provide  considerable  support 
for  the  form  proposed  in  Eq.  (5. 32). 

6.  Mass  and  Salt  Transfers  and  Halocline  Depths  in  an  Estuary. 

With  the  background  of  Section  5,  we  wiU  now  consider  the  effects  of 
finite  mixing  on  the  various  features  of  estuarine  circulations  and  distributions. 

We  have  in  mind  a body  of  water  like  the  Baltic  Sea,  although  the  model  has 
variable  parameters  which  permit  application  to  almost  any  estuary  whose 
mouth  has  a width  small  compared  to  the  general  horizontal  dimensions  of  the 
estuary.  A comparison  is  also  possible  with  laboratory  experiments  designed 
to  produce  estuarine  circulations. 

1 

t 
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We  will  not  discuss  at  any  laigth  the  practical  importance  of  studies  of 
estuarine  circulations  and  mass  distributions.  The  problem  of  pollution  of 
waters  near  the  surface  obviously  demands  attention,  but  another  problem  of 
strong  interest  for  the  Baltic  and  many  fjord-type  estuaries  of  Scandinavia  is 
the  stagnation  of  the  deep  water  (Fonselius,  1962,  1967,  1969;  Gade,  1970). 

This  is  associated  with  the  blocking  effects  of  the  sill  and  other  obstructions  to 
deep  water  motion  in  conjunction  with  the  stabilizing  effects  of  the  typical 
density  increase  with  depth  in  an  estuary  associated  with  temperature  and 
salinity  variations  in  the  vertical.  The  water  motions  are  very  slow  and 
therefore  even  exceedingly  weak  density  increases  with  depth  greatly  inhibit 
vertical  motions.  The  water  stagnation  implies  very  weak  or  zero  turbulence 
in  the  deep  waters  and  little  turbulent  mixing  with  the  fluid  above.  This  in  turn 
cuts  off  the  supply  of  oxygen  from  the  aerated  surface  water,  and  the  deep  layers 
may  ultimately  become  completely  exhausted  of  oxygen.  This  has  been  a pro- 
gressive development  in  the  deep  Baltic  over  the  past  75  years.  One  station  in 
the  central  Baltic  shows  an  oxygen  saturation  decrease  from  30%  to  near  0% 
during  this  period  at  a depth  of  160  m (Fonselius,  1969).  Fonselius  has  warned 
that  the  Baltic  deep  water  may  soon  become  devoid  of  life. 

The  complexities  of  real  estuaries  have  been  emphasized  by  many  authors, 
e.  g.  Pritchard  (1956),  Fonselius  (1969),  Gade  (1970)  and  Welander  (1974).  The 
Baltic,  for  example,  is  connected  to  its  source  of  salt  water,  the  Kattegat,  by  the 
Danish  Belts  and  the  Oresund.  This  in  itself  is  a great  complication  compared  to 
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only  the  Store  Belt  which,  in  fact,  is  mainly  responsible  for  the  fluxes  into  and 

out  of  the  Baltic.  We  take  the  minimum  width  W = 15  km  for  the  Belt  and  a mean 
depth  at  this  section  h = 18  m.  In  the  model,  the  estuary  waters  and  the  flow  at 

the  mouth  are  divided  into  two  layers  by  a thin  halocline.  This  is  not  a bad 
approximation  for  the  Baltic  which  has  a permanent  halocline  of  thickness  10  m 
at  a depth  D ^60  m,  although  there  certainly  are  vertical  and  horizontal  gradients 
of  salinity  above  and  below  the  halocline1 . For  example,  the  salinity  of  the  in- 
flowing water  in  the  Store  Belt  is  about  17.  5°/00  compared  to  ll°/00  in  the  Baltic  sea 
deep  water  although  the  model  ignores  this  and  assumes  no  mixing  as  the  ocean  water 

pours  into  the  basin.  The  salinity  of  the  upper  layer  in  the  Baltic  is  about  7°/00. 

The  model  assumes  that  the  inflowing  water  is  all  ocean  water;  by  way  of  com- 
parison, some  of  the  outflowing  Baltic  water  mixes  with  the  Skagerak  water  and 
recirculates  back  into  the  Baltic. 

The  water  balance  of  the  Baltic  region  involves  a close  equality  of  precipi- 
tation and  evaporation  so  that  the  fresh  water  influx  R is  equal  to  the  river  runoff 
which  we  take  to  be  R = 1. 49  x lO10  cm3 /sec.  The  outflowing  water  has  a flux  qx 
which  is  about  twice  the  magnitude  of  R (Brogmus,  1952). 

Estuarine  circulations  are  usually  very  small,  of  order  of  centimeters  per 
second  or  less,  and  this  suggests  a considerable  importance  of  friction.  This 
has  been  confirmed  by  many  investigations  (e.  g.  Pritchard,  1956)  and  is  treated 
in  Section  7.  At  a narrow  mouth,  however,  speeds  increase  very  considerably, 
inertial,  pressure  and  gravity  forces  become  comparable  (Long,  1975b)  and 
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friction  becomes  relatively  much  less  important.  In  the  present  Section  we  confine 
attention  to  flows  near  the  mouth  and  we  neglect  friction  entirely  in  this  region. 

As  we  discussed  in  Section  f,  we  adopt  the  mixing  law,  Eq.  (5. 32) 


where  we  estimate  K„  ~ 0.093  from  experiments  (Long,  1975b)  and  ctu  is  the  rms 
horizontal  turbulent  velocity  near  the  interface.  It  is  instructive  to  compute  the 
entrainment  velocity  for  two  choices  of  the  quantities  in  Eq.  (6. 1).  In  the  Baltic, 
for  example,  if  we  use  the  salinity  difference  of  4°/0o*  we  &et  &>  3.  2 cm/sec2 

The  turbulent  intensity  is  very  uncertain  and  u„  is  sensitive  to  its  value,  but  if 
we  use  an  estimate  ou  = 1. 4 cm/sec  suggested  by  a discussion  below,  and 
D = 60  m,  we  get  u,  =1.1  cm/day.  This  suggests  very  slow  adjustments  to 
changes  in  fresh  water  supply  and  turbulent  intensities.  In  estuaries  in  which  tidal 
effects  are  large  a„  may  be  much  larger.  If  we  use  oru  = 10  cm/sec,  D = 15  m, 

Cb  = 10  cm/sec2,  more  appropriate,  say,  to  estuaries  in  British  Columbia 
(Pi.kard  & Rodgers,  1959),  we  obtain  u,  = 5 m/day  and  adjustments  to  changes 
will  be  very  rapid. 

The  entrainment  velocity  in  the  estuary  may  be  written  u,  = q0/A  (Fig.  6. 1),  or  A 
is  the  area  of  the  halocline  in  the  estuary.  Using  Eqs.  (4. 14)  and  (4. 15),  we  have 


(6. 2) 

(6.  3) 
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We  may  now  relate  this  to  the  discussion  in  Section  4.  The  results  in 
Eq.  (4.  28)  apply  to  the  case  of  zero  mixing,  i.  e. , M = °°.  The  results  in 
Eq.  (4.  29)  imply  crtt  - « or  M = 0.  Obviously  then,  0 = AH0/h^  .equal  to 


8 = f(l-7) ) , M = <*> 


B = 1(1-  |-71),M  = 0 (6.6) 

in  the  two  extreme  cases  .will  be  some  function  of  7)  and  M for  finite  mixing. 

I don’t  see  any  way  to  find  this  relationship  rigorously  but  a simple  form  which 
satisfies  the  two  extremes  is 


3 = ! ri  - ■ 


3[l+M'f(7D] 


where  s is  a constant  and  f(71)  is  arbitrary.  There  are  additional  requirements 
which  must  be  satisfied.  Thus  as  cr„  - 0,  M - °°,  it  follows  from  Eq.  (3.  5)  that 


Eqs.  (4. 13)  and  (4.  25)  lead  to 


so  that  (6.  7)  and  (6.  9)  lead  to 


r-- 1(1-71) 


=3-371-20 


Using 
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(6. 10) 
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Eq.  (6. 10)  becomes 


s^i+.jgh*w; 


1 _6 


L1+l‘ 


) 


*01) 


(6. 12) 


D &T\3(tt>fW*  ' AK„< 

As  cru  - 0 all  other  quantities  remain  finite  and  we  see  that  s =2,  Ab  - b0.  Using 
(6.  8) , we  obtain 


f(Tl)  = ftfM)2 


We  get 


p = f1  _ ‘n[i6+27Maq-,n  )V  1 ( 

3[4+9M2(l-Tl  )2Tf  ] 


(6.13) 


(6. 14) 


A combination  of  (4.  l6),  (6.  2),  the  definition  of  F^,  (6.  9)  and  (4. 17)  leads  to 


Ab 


Q7  = 


\b. 


loi. 


; (1-7)  )3[4+9M2(l-TD2Tl4-47l] 


[4  + 9M2  (1 — T)  )2  if  J 


(6.15) 


D*  is  given  by 


V1  bn  - 


4TI4 


'0/'  (1-T)  )3  [4+ 9Ma(l-T))a  T)4-4T]] 


D*  = 


(6. 16) 


(6.17) 


Equations  (6. 14)-(6. 17)  determine  all  unknowns,  given  the  fres  water  flux 
ratio  Qf  and  the  mixing  number  M. 

The  form  of  the  solution  is  rather  complicated  but  a simple  example  of  a 
solution  of  Eqs.  (6. 14)-(6. 17)  is  an  estuary  with  an  infinitely  deep  sill,  i.  e. , 


h - ® . The  solution  then  depends  on  a single  non-dimensional  number 


fi.  # 

M'  = — — bQ-  z 
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Fig.  6.  2 Interface  and  halocline  depths  and  salinity  differences 
for  an  infinitely  deep  estuary.  D*  and  are  defined 
in  Eq.  (6.  l9).  is  the  same  as  £b/bg. 
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Using  the  definitions 


1 _ 2 1 _ 2 
v3u/3  u r'S 


, _ Dbn  W _ hi  bp  W 

if  ~ f “1  - ti 


R‘ 


R' 


where  hx  is  the  depth  of  upper  layer  at  the  constriction.  We  get 

1_ 

M'(l- 


n,  _ _ 1 v,  2 

~ Ah  * “1  (1—  ST  ) ~ 

M'n-^S)  D°  b0 


3M'hV 


IUI^ 


(6.19) 


(6.  20) 


Results  are  shown  in  Fig.  (6.  2).  The  depth  of  the  interface  at  the  mouth  hx 
and  the  halocline  depth  D increase  monotonically  as  M'  decreases  (mixing  in- 
creases). D is  about  twice  as  large  as  hx  for  larger  M'  over  most  of  the  range 
of  M'.  The  salinity  of  the  outflowing  fluid  increases  monotonically  as  the  mixing 
increases  and  the  fresh-water  flux  decreases.  The  flux  of  the  outflowing  fluid 
changes  very  slowly  with  M'  over  most  of  the  range,  qx/R  = (£b/b0)-1  increasing 
from  1.  0 to  1.  5 as  M' decreases  from  infinity  to  M'  ~ 1.  Below  M'  =1,  q X/R 
increases  more  rapidly  but  only  begins  to  exceed  2.  0 when  M'  drops  below  M'  = 0. 1. 

When  the  sill  depth  is  finite,  Eqs.  (6. 14) -(6. 17)  reveal  a common  behavior  for 
all  values  of  M,  namely  that  the  halocline  is  very  deep  for  both  small  and  large 
values  of  the  fresh-water  flux  R,  with  a minimum  at  a value  of  Qf , determined 
by  M.  The  quantity  Qf  B decreases  with  M;  for  example  Qf  „ =0.  244  at  small  M, 

QfB  = 0. 08  at  M = 100  and  QfB  = 0.  03  at  M = 1,000.  The  behavior  of  D*  as  a 
function  of  Qf  is  shown  in  Fig.  (6. 3)  and  (6. 4)  for  smaller  and  larger  values  of 
M,  respectively.  When  M exceeds  100  or  so  there  is  little  change  in  the  behavior 
over  most  of  the  range  of  Qf  except  near  Qf  B.  There  the  halocline  depth  continues 
to  decrease  strongly  as  M gets  still  larger. 


Fig.  (6.  5)  shows  an  experiment  by  Welander  (1974)  similar  to  the  model  of 
this  paper.  Data  are  lacking  to  permit  a careful  comparison, but  M = 8 seems  a 
good  choice.  Constants  were  chosen  to  force  agreement  near  the  minimum  point. 


0 0.1  0.2  0.1  0.4  OS  O.f  O.r  04  O.S  10 

0| 

Fig.  6.  3 Halocline  depth  as  function  of  fresh-water 
influx  for  smaller  values  of  M. 
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The  agreement  is  quite  good  except  for  small  values  of  R where  entrainment 
velocities  are  very  small  and  the  experiments  may  not  have  been  in  a steady 
state.  An  additional  factor  is  the  decrease  of  ctu  near  the  interface  as  D in- 
creases and  the  interface  gets  further  and  further  away  from  the  source  of 
energy  (stirrers).  This  effect  is  not  included  in  the  theory  and  would  tend 
to  reduce  D.  The  tendency  for  a minimum  depth  of  the  halocline  has  been 
observed  by  Tully  (1949)  in  Alberni  Inlet  as  we  see  in  Fig.  (6.  6).  The  theo- 
retical curve  corresponds  to  M = 0. 

The  quantity  0 varies  with  Qf  and  M as  seen  in  Fig.  (6.  7).  It  is  interesting 
to  note  that  gAS*  AH0/1T  varies  very  little  with  M.  Indeed,  one  would  anticipate 
that  this  quantity  would  vary  more  sensitively  with  Qf. 

The  influx  of  ocean  water  is  negligible  for  M > 100  or  so.  As  mixing 
increases,  Q0,  defined  by 


bo  h 2W 

increases.  Q0  rises  with  increase  of  Qf , reaches  a naximum  and  then  de- 
creases. Examples  are  shown  in  Fig.  (8.  8).  The  curves  for  M < 3 or  so  do 
not  differ  very  much  from  the  curve  shown  for  M = 0. 10.  This  typical  behavior 
of  Q0  appears  in  other  theories  of  es+ua-ies,  for  example  Kullenberg  (1955). 

The  salinity  of  the  outflowing  fluid  increases  as  the  fresh  water  influx 
decreases.  This  was  also  pointed  out  by  Kullenberg.  Fig.  (6. 9)  shows  the 
curves  for  two  values  of  M.  Decreases  in  mixing  intensities  accompany  de- 
creases in  salinity,  as  expected. 

An  estimation  of  M for  various  applications  is  important  but  M is 
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extremely  sensitive  to  ou  and  au  is  difficult  to  estimate.  We  may  m'dce  the 
following  computations  for  the  Baltic  Sea:  Fonselius  (1969)  gives  the  ratio 
qj,  /R  "=  2.  We  may  compute  Qf  from  Eq.  (4. 10)  with  the  constant  equal  to 
-gAH0/6b  using  the  data  of  Fonselius  of  R = 1.49  x l0locm3/sec  and  width 
W = 15  km,  depth  h = 18  m for  the  Store  Belt.  The  appropriate  value  of  b0 
corresponds  to  a salinity  in  the  lower  layer  of  17.  Ef/ '00  or  b0  = 14cm/ sec3. 

This  leads  to  Qf  = 0.  0348.  Finally  D/h  for  the  Baltic  is  60/18  or  3.  33.  If 
we  use  the  last  two  figures,  we  obtain  close  agreement  for  a value  of  M = 12. 
The  computations  yield  Qf  = 0.  0346  and  D/h  = 3.  45.  The  value  of  q±  /R,  how- 
ever, is  3.  29  which  is  considerably  too  large.  The  value  of  M = 12  permits  us 
to  compute  cru  using  A = 3. 1 x 101  5cm2.  We  obtain  ou  =1.4  cm/sec  which  is, 
perhaps,  reasonable. 

One  may  speculate  regarding  the  changes  in  the  Baltic  over  the  past  75 
years  as  described  by  Fonselius.  The  stability  of  the  Baltic  has  apparently 
increased,  because,  although  the  salinity  of  the  upper  and  lower  layers  have 
both  increased,  the  salinity  of  the  lower  layer  has  increased  somewhat  more. 
This  stability  increase  may  be  related  to  the  reduction  of  oxygen  in  the  deep 
Baltic.  R has  decreased  about  15%  over  this  period  and  Fonselius  ascribes  the 
increase  of  salinity  in  the  Baltic  as  a whole  to  this  decrease  of  R.  A further 
change  has  been  a decrease  in  the  depth  of  the  halocline  from  80  m to  60  m over 
the  period.  According  to  present  theory,  for  the  value  of  M = 12  and  the  other 
conditions  of  the  Baltic,  a decrease  of  R with  all  other  factors  held  constant 
should  have  resulted  in  an  increase  of  D.  However,  as  we  have  seen,  D is  most 
sensitive  to  (Ju.  Since  precipitation  and  runoff  have  decreased,  this  implies  a 
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Fig.  6.  8 Variations  of  ocean  water  flux  with  fresh  water  influx. 
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decrease  of  storminess  over  the  Baltic  as  asserted,  in  fact,  by  Fonselius,  and 
a decrease  of  <ju  . This  could  have  been  overriding  and  serves  to  explain  the  de- 
crease of  D.  Notice  that  only  long  period  changes  in  au  are  involved  because  we 
have  decided  that  the  response  of  the  halocline  depth  to  changes  is  very  slow. 

7.  Circulations  and  Density  Distributions  in  a 


In  Section  6 we  considered  an  estuary  of  a rather  general  type  except  that 
the  connection  to  the  open  sea  has  such  a character  that  (1)  the  flow  in  this  portion 
of  the  estuary  (mouth)  is  perfect  fluid  flow  and  (2)  conditions  become  critical  at 
one  section  of  the  mouth.  Our  investigation  derived  an  expression  for  the  depth 
of  the  halocline  in  the  estuary  proper  but  other  than  this  nothing  was  found  out 
about  the  circulations  and  salt  distributions  in  the  estuary  except  near  the  mouth. 
We  now  investigate  (Long,  1975b)  conditions  throughout  an  estuary  without  the  two 
assumptions  mentioned  above.  For  simplicity  we  assume  that  the  estuary  is 
infinitely  deep  everywhere.  As  it  develops,  this  means  that  the  halocline  depth 
is  small  compared  to  the  depth  of  the  estuary.  This  is  not  a bad  assumption  in 
many  cases,  for  example  in  fjords  in  British  Columbia  and  Alaska,  but  the  model 
does  not  include  the  situation  in  some  Scandanavian  fjords  in  which  there  is  a sill 
whose  depth  is  of  the  order  of  the  depth  of  the  halocline. 

The  model  (Fig.  7. 1)  has  a well-mixed  upper  layer  and  a deep,  non-turbulent 
uniform  lower  layer  of  the  density  of  sea  water.  The  width  of  the  estuary  is  W(x) 
and  the  sides  are  vertical.  The  free  surface  is  given  by  z = H and  the  interface  by 
z = D0(x).  The  upper  layer  has  a mean  depth  D(x).  The  variation  of  H along  the 
estuary  provides  an  important  component  of  the  pressure  gradient  force  but  incre- 
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merrts  in  H are  of  the  order  of  10  centimeters  over  tens  of  kilometers  and  so  may 
be  neglected  elsewhere  in  the  argument  and  H taken  to  be  a constant.  The  vertical 
variation  of  the  buoyancy  is  confined  to  a thin,  interfacial  zone  whose  thickness  is 
neglected,  but  the  mean  buoyancy  varies  continuously  with  x in  the  upper  layer  as  the 
water  flowing  out  to  sea  becomes  more  and  more  brackish.  In  general,  a symbol 
such  as  b denotes  an  ensemble  average  or  average  over  a long  time,  and  b denotes 
an  average  of  b over  the  vertical  cross-section  in  the  upper  fluid. 

The  buoyancy  difference  between  the  two  layers  is  Ab  and  varies  along  the 
channel  because  the  buoyancy  of  the  upper  layer  varies.  The  rms  turbulent 
velocity  in  the  upper  layer  is  ou  and  we  assume  that  ou  is  uniform  in  the  layer. 

To  be  definite  we  take  o„  to  be  the  rms  horizontal  velocity  component  along  the 
x-axis.  The  non-dimensional  quantity  Rig  = D£b/o^  has  the  form  of  a Richardson 
number.  If  we  take  as  typical  values,  D = 103  cm,  Ab  = 25  cm  sec  , <7U  =10  cm  see  , 
we  get  Ri2  = 250.  We  have  taken  a rather  high  value  for  ou  so  that  Ri2  will  usually 
be  greater  than  this.  It  is  reasonable  then  to  assume  that  the  estuary  is  strongly 
stratified  in  the  sense  that  Ris  » 1. 

Conservation  of  mass  for  a region  of  the  upper  layer  of  length  Ax,  width  W, 
and  thickness  D yields 

A(WDu ) = JJVn  da  (7.1) 

where  uis  the  mean  velocity  over  the  cross-section  of  the  upper  layer,  da  is  an 
element  of  area  of  the  interface,  and  the  integral  is  over  the  interfacial  boundary 
of  the  region.  The  normal  velocity  into  the  layer,  vn  , is  zero  when  there  is  no 
mixing.  We  have  already  assumed,  however,  that  turbulence  exists  in  the  upper 
layer,  and  we  may  identify  vn  with  the  entrainment  velocity  so  that  with  use  of 
Eq.  (5.32),  Eq.  (7.1)  becomes 
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^(WDu) 


(7.2) 


The  pressure  in  the  lower  layer  is 

? = g(H-z)  + b(H-D0)  + b0(D0-z)  (7.  3) 

Pr 

If  we  assume,  as  indicated  above,  that  the  lower  fluid  has  a depth  large  compared 
to  that  of  the  upper  fluid,  horizontal  accelerations  will  be  negligible  and  the  longi- 
tudinal pressure  gradient  is  zero  in  this  layer  (Gade,  1974).  Then 
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Let  us  now  integrate  the  horizontal  equation  of  motion, 


v.  uv 


pH  + m-zl  — I __§IU'U')  _ JL(U'V')  _ ® (u'W)  5. 

gH*+(Hz)dxJ  dx  dy  dz  { ,5) 


over  a region  of  the  upper  fluid  of  length  Ax,  width  W and  thickness  D.  We  get 
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(7.6) 


where  the  pressure  in  the  upper  layer  is  given  by  p = ft g(H-z)+  pfb(H-z)  and 
where  t,  is  very  nearly  equal  to  the  average  stress  at  the  interface.  We  assume1 
Tt  = Ku3  where  K is  the  drag  coefficient.  A term  involving  the  horizontal  rate  of 
change  of  tx  has  been  omitted  in  Eq.  (7.  6)  because  the  turbulence  is  assumed  to 


xThis  is  a standard  assumption  for  example  in  flow  in  pipes.  The  stress  is 
basically  equal  to  -TTwT  and  in  typical  shearing  flows  ou/u  and  uM/u  are  small  and 
vary  very  slowly  with  Reynolds  number  and  so  are  taken  to  be  constants.  The  stress 
then  is  proportional  to  a with  a small  coefficient  of  proportionality.  In  the  present 
case  it  is  likely  that  cru/l  and  %/xi  are  of  order  one  so  that  K is  of  order  one.  We 
assume  this  later. 
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be  uniform  in  the  upper  layer.  If  ut  is  the  average  horizontal  velocity  at  the 
bottom  of  the  upper  layer,  the  integral  on  the  ihs  of  Eq.  (7.  6)  may  be  written 

%^J,J’uudzdy  _U*Ul  (7‘7) 

The  ratio  of  ue  u5  to  the  term  Tt  in  Eq.  (7.  6)  is  of  order  u,  D/K„ , or  less,  where 
K„  is  the  eddy  viscosity  in  the  upper  fluid  near  the  interface.  We  estimate 
Kb  as  £cru  where  JL  is  the  eddy  size.  If  we  use  observations  of  flow  in  pipes,  we 
may  estimate  D/3  (Hinze,  1959).  Using 

~ _1_  of 

u«  = 10  DAb 

the  ratio  is  approximately  0. 3 cr„/Dib  which  we  have  assumed  to  be  small. 

Thus  the  second  term  in  Eq.  (7.  7)  may  be  neglected. 

Let  us  now  write 

JJuudzdy  = yuuDW  (7.  8) 

The  quantity  y depends  on  the  velocity  distribution  in  the  layer.  If  the  velocity 
is  uniform,  y = 1.  Other  physically  reasonable  assumptions  yield  values  greater 
than  1 but  less  than  2 or  so.  Our  analysis  does  not  yield  a precise  value  for  y, 
but  we  find  that  y occurs  in  combination  with  other  equally  uncertain  constants 
of  the  problem  so  that  it  seems  pointless  to  attempt  to  refine  this  portion  of 
the  argument. 


Eq.  (7.  6)  becomes 
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]_The  neglected  quantity  is  of  order  u'b'D.  The  ratio  of  this  to  ubD  is 

u'b'/ub  ~ub/b  since,  as  we  discuss  later,  ffu  ~u.  If  the  horizontal  eddy  dimension 
is  of  the  order  of  the  width,  ob  is  of  the  order  of  the  longitudinal  mean  buoyancy 
variation  over  a distance  W.  This  is  small  compared  to  b if  the  estuary  is  long 
compared  to  the  width. 
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7.  2 The  Solution  for  an  Estuary  of  Uniform  Width. 


Let  us  consider  the  solution  of  (7. 19)  when  X = 1 , i.  e. , when  the  estuary 
has  a uniform  width.  If  we  use  a new  dependent  variable  C = if  /Q3,  and  a new 
independent  variable,  inT!  , we  obtain 


3 T 3mIC  6m/ 


ip1 — ~ ^ =d£nt1 

+ + £ 


The  integral  of  this  is 


i+_JL_  \ 

3 3mK./ 


mK„  / 


. , _K_  j " 3 £nV1+  mK„+  2m/1  = m T'  + const 

mK,,  + 2m 


Using  new  constants, 


s — 1 + 


, C3  = 


l/2m 


(7.20) 


the  solution  is 


T)  = " rrr-  + const 

(1  +C3  Q ~T~ 


(7.21) 


Imposing  the  boundary  condition  that  7)  = at  the  head  of  the  estuary  where 


Q =1,  Eq.  (7.  21)  becomes 


H = i -£j ' i3  ItgjTg,  3 
TI  \n?J  si  +&r  J 


\ si+c^c  ) 

The  quantity  ? may  be  found  by  integrating  Eq.  (7. 18).  We  get 

? = P Q d"J 

Q* 

where  we  choose  Q*  to  be  the  flux  at  an  arbitrary  section  at  which  e = 0. 


(7.  22) 


(7.  23) 


r 


We  may  simplify  the  problem  considerably  by  the  new  definitions 
C =C3C.  TV  = CQJT1 , Q'  =Q;Q,  r = 2^nCQh'| 


(7.24) 


where  we  denote  by  the  value  of  Q'  at  the  head  of  the  estuary.  We  then  obtain 


-ta"  jii3 

* C = 


(1+  C')  3 ** 

Q 71!  t~*  ® **  1 77  6 ■- 1 

tUq\  Qi=- »+i 

q;  q (i+c3^3)  3 


(7.25) 


(7.26) 


where  Q*  = QiQ*.  Notice  that  1 < s < 2 tending  to  s = 1 as  cru  decreases  to  zero 
and  to  s = 2 when  <ru  is  very  large  or,  on  the  other  hand,  as  friction  becomes  more 
or  less  important,  respectively. 

We  may  solve  Eqs.  (7.  25)  for  Q\  We  get 


(i+  o 


Differentiating  Eq.  (7.27),  we  obtain 


dQ-  I \rr<W^ 


(7.27) 


(7.  28) 


so  that  Q'  has  a maximum,  q;,  when 

r ' - S~1 

to  - 2 

The  corresponding  maximum  value  of  Q is  Qe.  Notice,  also,  that 


(7.  29) 


C3dQ'  =d?' 


(7.  30) 


so  that  £'  and  Q'  increase  or  decrease  together.  This  simply  means  that  the  flux 
in  the  upper  level  increases  seaward,  as  is  obvious  for  physical  considerations. 
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Since  F is  a maximum  when  Q is  a maximum,  the  flux  may  increase  to  the  maximum 
as  we  move  seaward  but  if  this  happens  the  estuary  of  constant  width  must  then  end. 

It  is  enlightening  to  compute  the  Froude  number  F at  Cj.  A convenient 
definition  is 


F8 


D& 


(7.31) 


This  yields 


(7.  32) 


Using  (7.  24)  and  (7.  29),  we  get  F^  = 1,  so  that  the  maximum  flux  Qe  corresponds 

to  a point  of  critical  flow  with  supercritical  flow  for  £'  < £c'  and  subcritical  flow 

for  C > Cc  • We  may  compare  this  to  a theory  by  Stommel  (1951)  which  led  to  the 

differential  equation  in  Eq.  (7. 19)  with  X = 1 and  with  the  frictional  term  missing.  He 
obtained 


U dQ  1-F  8 
Q dl)  “ 2(^-F2) 


(7.  33) 


Stommel' s theory  was  incomplete  and  he  could  not  solve  for  fluxes  and  interface 
depths  as  functions  of  distance  along  the  estuary. 

From  Eq.  (7.  25)  we  get 


•—3 

dJL  = 3 (s-C) 

dc'  6+4 

(1+0  3 


(7.  34) 


Let  us  now  consider  the  general  properties  of  the  solution.  If  we  increase  £' 
from  £e'  = (s-l)/2  (corresponding  to  subcritical  flow),  s-C'  is  positive  at  first 
and  7i'  increases.  It  reaches  a maximum  at  £'  = s and  then  decreases  as  £'  gets 
larger.  Since  Q'  is  decreasing,  Eq.  (7.30)  shows  that  F'  is  decreasing  so  that 
we  are  moving  toward  the  head  of  the  estuary.  On  the  other  hand,  if  we  decrease 
C'  from^  (corresponding  to  supercritical  flow),  we  find  that  7)  ' decreases  mono- 
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tonic  ally.  Again  %'  decreases  so  that  we  are  again  moving  toward  the  head  of  the 
estuary.  Thus,  two  distinct  solutions  are  possible,  but  the  second  corresponds  to 
high  velocities  in  excess  of  a meter  per  second  or  more  at  the  head  of  the  estuary 
and  elsewhere  and  seems  unlikely  to  occur.  It  also  involves  a decrease  of  velocity 
and  increase  of  elevation  of  the  free  surface  toward  the  mouth  and  both  are  contrary 
to  observations  (Gade,  1970).  Consequently,  we  adopt  the  solution,  C'  •>  Gc»  f°r 
subcritical  flow. 

The  remaining  problem  concerns  conditions  at  the  mouth  of  the  estuary. 

Let  us  assume  that  the  estuary  has  a uniform  width  from  the  mouth  to  the  head 
but  at  the  mouth  the  width  suddenly  increases  rapidly.  The  arguments  leading 
to  Eqs.  (7. 14)— (7. 16)  permit  a variable  width.  Accordingly,  let  us  derive  two 
equations  from  Eqs.  (7. 15)  and  (7. 16)  yielding  the  rate  of  change  of  H and  of 
F8  = mQ3/Tp>2  along  the  channel.  We  get 


dTl_  _1_  2K„m\  1 KF8  -a  ]]F8d\ 

df  - l-F2  ^ ” V A~  ->TY,V  \ ~ * 1 Ti  rlF  I 


W"  2mK„X  r X d§J 


(7.  35) 

dF2  3F2K„mXr^.  1 . KF2  ~ F2(F2  + 2)  dX  „ 

df  = # f)'Ln2+  mK^  ~ '\(1-FT  ^ (7-36) 

Let  us  now  investigate  all  possibilities,  shown  in  Fig.  7.2,  for  the  location  of 

the  mouth  of  the  estuary. 

(a)  The  estuary  ends  at  or  before  the  section  of  maximum  depth  of  the 


upper  layer  (point  I or  J). 

In  this  case  1 - F2  > 0 and  the  quantity  in  curly  brackets  in  Eq.  (7.  35)  is 
positive  to  the  left  of  I or  J so  that  T)  is  increasing  with  |.  Eq.  (7. 36)  reveals 
that  F2  is  increasing  but  is  considerably  less  than  one.  Just  past  point  I or  J, 
where  X is  still  close  to  one  but  dx/d|  is  very  large,  the  last  term  in  curly  brackets 
in  Eq.  (7.35)  will  quickly  dominate  and  11  will  increase  even  more  rapidly.  At 


A 


f 
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the  same  time,  the  last  negative  term  in  Eq.  (7. 36)  will  dominate  and  F2  will 
start  decreasing.  If  we  now  regard  as  physically  necessary  that  T]  should  de- 
crease past  the  mouth  as  the  mass  of  brackish  water  spreads  in  all  directions, 
this  case  is  impossible. 

(b)  The  estuary  ends  after  the  section  of  maximum  depth  of  the  upper 
layer  but  not  close  to  the  section  at  which  F2  = 1. 

In  this  case,  T|  is  decreasing  just  to  the  left  of  point  K and  F2  is  somewhat 
smaller  than  one  and  increasing.  Just  past  K,  however,  the  last  term  in  curly 
brackets  in  Eq.  (7.  35)  quickly  dominates  and  71  will  reach  a minimum  and  then 
increase.  Also,  the  last  negative  term  in  Eq.  (7. 36)  begins  to  dominate  and 
F2  begins  to  decrease.  Again  it  is  not  possible  for  71  to  decrease  past  the  mouth 
and  this  case  is  also  impossible. 

(c)  The  estuary  ends  just  before  the  section  at  which  Fa  = 1. 

In  this  case,  71  is  decreasing  to  the  left  of  point  L and  F3  is  very  close 
to  one.  Then,  if  dTl/d?  remains  finite,  the  quantity  in  curly  brackets  in  Eq. 

(7.  35)  falls  to  zero  at  a section  just  past  point  L,  F2  passes  through  the  critical 
value  and  1 - F2  becomes  negative.  After  this,  the  last  term  in  curly  brackets 
in  Eq.  (7,  35)  dominates  and  7]  continues  to  decrease.  Eq.  (7.  36)  shows  that  F2 
also  continues  to  increase.  This  behavior  is  entirely  reasonable  and  we, 
therefore,  accept  the  point  L as  the  end  of  the  estuary.  It  is  convenient  to  take 
this  section  as  the  origin  of  our  coordinate  system  so  that  § = §'  =0  there  and 
in  Eqs.  (7.  23)-(7. 26)  Q+  and  may  be  replaced  by  Q0  and  . 

7.  3 Discussion. 

We  see  from  the  theory  that  7)’  = g(?',s)  or,  from  Eq.  (7.24),  C71Q{,  = g(§',s). 


i 


r 


k 


namely  the  inner  Oslofjord  (Gade,  1970)  and  the  Knight  Inlet  (Pickard  and  Rodgers, 
1959),  for  which  we  use  qf  = 60  cm3/sec  and  2.  3 103  cm2/sec  respectively. 

It  is  important,  of  course,  to  estimate  m in  Eqs.  (7. 14),  but  this  quantity  is 
extremely  sensitive  to  the  value  of  cr„.  In  either  case,  however,  and  in  most 
estuaries,  it  seems  likely  that  m is  considerably  less1  than  one.  The  quantity 


has  been  estimated  as  K„  = 0. 1.  The  constant  K is  uncertain  but,  as  we  have 
indicated,  we  will  use  K = 1.  The  estimates  we  have  made  so  far  mean  that  s 
is  close  to  but  a little  larger  than  1.  The  computations  from  the  theory  reveal 
little  change  in  Q0  as  s -*1.  At  mKnCLo^/qfy  = 3,  for  example,  Qc  increases 
by  about  10%  as  s decreases  from  the  rather  large  value  of  1. 2 to  s = 1.  0036. 
According  to  Fig.  7.3,  we  may  write,  approximately, 

1 In  Knight  Inlet,  for  example,  we  may  take  y = 1. 3,  =5  cm/sec, 

b0  = 25  cm/sec3,  qf  =2.3  103  cm3 /sec.  Then  m = 2.  8 10  . In  the  very  different 

Oslofjord,  we  may  take  y = 1.3,  cru  =1  cm/sec,  b0  = 5 cm/sec-3,  qf  = 60  cm3/sec. 
Then  m =4.3  10  . 


mKn  CLou 

Yq, 

Relationships  between  the  halocline  depth  at  the  hea 
of  the  estuary  and  the  estuary  length. 


rnKnCLffu 

Yq, 

Fig.  7.  3 Relationships  between  the  mass  flux  in  the  upper  laye: 
and  the  estuary  length 


Fig.  7.4 
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Qe  = 1+0.9 


(7.  37) 


Our  estimates  lead  to  C =0. 37  and,  with  y =1.3,  Eq.  (7. 37)  may  be  written 


Qe  *=  1 +0.  033- 


(7.38) 


In  the  Oslofjord,  with  L ”=  2.  5 106  cm,  b0  = 5 cm/sec3  and  Q0  *=  3. 2 (Gade,  1970),  we 
estimate  c ru  = 0.  83  cm/sec.  In  the  Knight  Inlet,  Qc  is  quite  uncertain,  but  the  data 
suggests  a value  of  perhaps  4.  Using  L ~l. 1 107  cm,  b0  ~ 25  cm/secs,  we 
obtain  ctu  = 5.  7 cm/sec. 

Computations  for  various  values  of  s near  1 reveal  little  variation  of  CT)  „ 
with  s for  all  but  very  short  estuaries.  A rough  relationship  is 


or,  approximately. 


CT)h  = 1 + 0.25 


D„  ~ 2.  7 |f  + 0.  025  L 


(7.  39) 


(7.40) 


Computations  in  Fig.  7. 3 and  7.  4 yield  CT^  =1.7  and  1.  9,  respectively,  for  the 
Oslofjord  and  the  Knight  Inlet.  We  get  approximately  Dh  = 3. 3m  for  the  Oslofjord 
and  21  m for  the  Knight  Inlet.  The  first  is  a considerable  underestimate  for  the 
Oslofjord  but  the  second  is  close  to  observations  in  the  Knight  Inlet.  The  com- 
parisons suggest  that  the  model  of  this  paper  may  contain  the  basic  physical 
mechanisms  of  the  Knight  Inlet  but  probably  differs  fundamentally  from  the  Oslo- 
fjord. This  is  not  surprising.  The  former  has  a geometry  similar  to  the  model 
whereas  the  latter  is  very  different.  For  example,  the  Oslofjord  has  a shallow  sill 
depth  which  forms  a considerable  barrier  for  the  influx  of  salt  water,  whereas  the 
sill  depths  of  the  Knight  Inlet  are  well  below  the  halocline.  There  is,  moreover, 
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the  observed  basic  difference  in  the  horizontal  density  variation.  This  is  similar 
to  the  model  in  the  Knight  Inlet  but  is  virtually  absent  in  the  Oslofjord,  so  that  the  • ] 

basic  driving  mechanism  of  the  present  model  is  absent  in  the  latter.  As  we  dis- 
cuss in  Section  8,  the  Oslofjord  appears  to  be  influenced  primarily  by  the  density 
distribution  in  the  waters  outside  of  the  fjord. 

Notice  that  if  we  consider  Dh  as  a function  of  qf , Eq.  (7. 40)  shows  that  the 
halocline  depth  is  a minimum  for  a certain  value  of  qf  and  large  for  both  small 
and  large  values  of  the  fresh-water  discharge.  This  is  the  same  as  the  behavior 
discussed  in  Section  6. 

An  interesting  feature  of  the  theory  is  the  slope  of  the  free  surface.  As 
we  have  noted,  the  surface  slopes  downward  toward  the  mouth  in  subcritical 
flow  and  Eq.  (7.4)  yields  for  the  total  drop  AH, 

AH  = -^-AC^  (7.41) 

A ^ 

For  the  Knight  Inlet,  for  example,  AC'3  = 2,  and  we  get  Afl  = 55  cm.  This  seems 
rather  large  but  values  of  10-15  cm  in  shorter  Norwegian  fjords  have  been  found 
by  Gade  (personal  communication).  Notice  that  the  commonly  used  argument  that 
the  free  surface  must  slope  downward  tow  rd  the  mouth  (Gade,  1974)  implicitly 
assumes  subcritical  flow.  We  have  not  discussed  the  supercritical  case  in  which 
the  free  surface  slopes  upward  toward  the  sea,  but  it  is  possible  that  some 
estuaries  have  this  character. 

The  calculated  drop  of  surface  level  for  the  Knight  Inlet  represents  a 
potential  energy  far  in  excess  of  the  observed  or  theoretical  kinetic  energies 
and  we  conclude  that  friction  dominates  flow.  This  is  probably  true  of  most 
fjord-type  estuaries.  We  also  deduce  the  dominate  effect  of  friction  from  the 
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fact  that  s tends  to  be  very  close  to  1.  As  we  see  in  Eq.  (7. 20)  this  corresponds 
to  large  values  of  the  drag  coefficient  K.  Indeed,  the  ratio  of  inertial  forces  to 
pressure  or  frictional  forces  is  proportional  to  s-1. 


7.  3 Estuary  of  Arbitrary  Width  - Integration  Procedure. 

We  may  outline  the  procedure  for  finding  the  solution  for  an  estuary  of 

arbitrary  shape.  We  first  notice  that  Eqs.  (7.  35)  and  (7.  36)  are  invariant  for 

* * $ 

transformations  T1  = orT)  , C = «?  . Using  X = f(§  ),  we  may  write 

= -^“-^-^(J-F2)-  -fFs+T1*F2  — * (7.4 

a V Y d5 


Tl*(l-F2)f  jjp  = 3^-m  F2 f2  (F2  + a')+  — ■ 


-Fs(Fa+  2)  1)*  (7.43) 


Let  us  investigate  the  solution  in  the  region  just  past  the  mouth  (F*=  0)  where,  by 
assumption,  the  width  of  the  estuary  increases  linearly.  At  5*  = a,  the  flow 
becomes  critical.  Then  if  £*  = 5*-  a and  \ is  the  value  of  X at  f*  = 0, 


= Tic*  + bi  C+  •••.  F2  =l+aiC*+...,  f = X.  + Cia+CiC*  (7.44) 


Substituting  (7. 44)  into  (7. 42)  and  (7.  43)  and  equating  coefficients  of  £ * we  obtain 
to  zero  order  two  identical  equations.  This  may  be  used  to  solve  for  X,+  cx  a: 


JS+.'C  frhK.mil? 


X,+  Ci.a  = 


V V 


(7.45) 


3K„  m/y 

The  other  two  equations  obtained  from  equating  coefficients  of  r*  may  be  written 


Aibi  + BiTl?  = Ej.  (7.46) 

Aabi  + BaT^sEa  (7.47) 


I 
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where 

A;  = ax  (X„+c1a)+c1 , Bx  =c1ai,  As  = -3cx , B2  = a12(X,+c1a)-4c1a1  (7.48) 
Ex  = + -^-^|c1(X.+c1a)  + a^X.+  ^a)2  + *[ax  (X.+cxa)+Cj.  (7.49) 

E2  =-^—  |a1(XB+cxa)2+3c1(>B  + c1a)j- ^ 2ax(Xll+c1a)+cx  ; (7.50) 

Thus 

b,  = (E1B2-EsB1)/(A1B2-AsB1)  (7.51) 

Tl*  = (Ax  E2-A2Ex  )/(A1B2-A2B1 ) (7.  52) 


The  integration  may  then  proceed  as  follows.  We  specify  X,  and  cx  and  the 
constants  Kn , m,  y.  We  also  specify  ax.  This  permits  calculation  of  bx  and  7|* 
from  (7.  51)  and  (7.  52),  and  X.+  ^a  from  ;7.45).  We  may  then  compute  7]*  and 
from  (7. 44)  for  say  F*  = 0,  and  F*  = -Ac*.  We  may  then  use  the  finite  difference 
forms  of  (7.42)  and  (7.43)  to  compute  F2  and  7]*  at  subsequent  grid  points  until 
X = 1.  This  corresponds  to  the  head  of  the  estuary,  where  r = Fh.  Q = 1, 

T),  =71*0.  Then  we  have 


3 

O 


m 


F2^*3 


(7.  53) 


and  we  have  a solution  for  an  estuary  of  width. 


X 


J 


(7.  54) 


8.  Three-Layer  Circulations  in  Estuaries  and  Harbors. 

A number  of  writers  beginning,  apparently,  with  Hachey  (1934)  have  remarked 
that  if  a density  variation  with  depth  exists  in  the  water  outside  of  an  embayment 
and  if  there  is  a source  of  mixing  in  the  embayment  from  tides  and  wind,  there  will 
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fact  that  s tends  to  be  very  close  to  1.  As  we  see  in  Eq.  (7. 20)  this  corresponds 
to  large  values  of  the  drag  coefficient  K.  Indeed,  the  ratio  of  inertial  forces  to 
pressure  or  frictional  forces  is  proportional  to  s-1. 


7.  3 Estuary  of  Arbitrary  Width  - Integration  Procedure. 

We  may  outline  the  procedure  for  finding  the  solution  for  an  estuary  of 
arbitrary  shape.  We  first  notice  that  Eqs.  (7.  35)  and  (7.  36)  are  invariant  for 

afc  afe 

transformations  Tl  = cvTi  • 5 = ckF  . Using  X = f(5  ),  we  may  write 


f(l-F2)~J*-  = -^^-^(J-F2)-  — fF2  + 71*  F2  — * 
a*1  Y Y de 


(7.42) 


*(1-F2)f  — * = — .a11?  F2 f2 (F5  + a)+  — f -Fa(Fa-f2)Tl* 


dF 


df 

dF* 


(7.43) 


Let  us  investigate  the  solution  in  the  region  just  past  the  mouth  (F*=  0)  where,  by 
assumption,  the  width  of  the  estuary  increases  linearly.  At  = a,  the  flow 
becomes  critical.  Then  if  £*  = 5*-  a and  \ is  the  value  of  X at  f*  = 0, 


Tl*  = TVs*  + biC^  ....  F2  a 1+ax  C*+. . . , { = X.  + Cia+CiC*  (7.44) 

Substituting  (7.44)  into  (7.42)  and  (7.43)  and  equating  coefficients  of  £*  we  obtain 
to  zero  order  two  identical  equations.  This  may  be  used  to  solve  for  X.+  Cj^a: 

-K  +/*£+ 

X.+  e^  = -- ‘1—  t * (7.  45) 

3K„  m/y 

The  other  two  equations  obtained  from  equating  coefficients  of  r*  may  be  written 

Aibi  + BiTl*  = Ex  (7.46) 

Asbi+B^  = EB 


(7.47) 
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I 


be  a tendency  for  a density  distribution  as  shown  in  Fig.  8. 1 and  a resulting 
three-layer  circulation  from  the  density-pressure  effect  even  in  the  absence 
of  a fresh-water  influx.  Stroup,  Pritchard  and  Carpenter  (1961)  have  calculated 
that  this  type  of  circulation  is  dominant  as  a flushing  mechanism  for  Baltimore 
Harbor  in  the  Chesapeake  Bay.  Here  the  less  salty  water  in  the  upper  layer 
outside  of  the  harbor  originates  from  the  fresh  water  discharge  of  the 
Susquehanna  River  at  the  head  of  the  Chesapeake  Bay.  It  has  also  been 
suggested  that  the  outside  density  distribution  is  important  for  the  circulations 
in  the  Oslofjord  in  Norway  (Gade,  1970)  and  in  the  Gullmarfjord  on  the  west 
coast  of  Sweden  (Rydberg,  1975).  Rydberg  observed  inflow  in  the  layer  below 
the  halocline  at  the  sill  and  inferred  from  considerations  of  conservation  of 
salt  that  there  must  be  an  influx  of  lower  salinity  water  somewhere  in  the  upper 
level.  Conservation  of  mass  requires  an  efflux  of  water  at  some  intermediate 
level.  The  lower  salinity  water  outside  of  the  Oslofjord  and  the  Gullmarfjord 
originates  as  the  brackish  water  flowing  out  of  the  Danish  sounds  from  the 
Baltic  Sea.  The  three-layer  type  of  circulation  has  also  been  discussed  by 
Hansen  and  Rattray  (1972). 

A simple  laboratory  model  of  a three-layer  circulation  was  constructed 
by  Hachey  (1934).  It  is  shown  schematically  in  Fig.  8.  2.  A long  channel  con- 
tains, initially,  two  layers  of  water  of  different  density.  At  one  end  of  the  channel 
a cylinder  at  the  interface  is  rotated  as  shown  to  provide  a source  of  mixing. 

The  resulting  pressure-density  distribution  produces  the  three-layer  circulation 


in  the  figure. 

The  simplest  model  of  the  basic  phenomenon  of  this  section  is  one  in  which 


producing  three-layer  circulation 
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there  is  an  ocean,  infinite  laterally,  composed  of  two  layers  with  buoyancy 
difference  £b0  and  with  an  upper  layer  of  depth  h.  These  are  the  conditions 

outside  of  an  embayment  of  arbitrary  geometry  specified  by  lengths  1^  , L2 

The  fresh  water  influx  is  zero.  There  is  some  source  of  mixing  in  the  estuary. 
The  nature  of  the  mixing  is  not  important  for  the  ultimate  argument,  I think,  but 
for  simplicity  let  us  consider  a mechanism  similar  to  the  one  we  use  in  the 
experiment  discussed  at  the  end  of  this  section,  namely  a rotating  propellor 
located  at  depth  h of  angular  speed  oj  and  geometry  specified  by  the  lengths 
ai  9 a2 , • • . . This  mixing  will  result  in  an  efflux  in  the  vicinity  of  the  depth  h or 
magnitude  qt.  Dimensional  analysis  yields 


-5--^  1 = f . ti.  |l2 Sl_  S2 


h7(Ab0)* 


Ab0  ’ h ’ h 


h ’ h *•* 


(8.1) 


Now  let  u increase  from  a low  value,  qj  increases  with  w.  At  some  value  of  w 
the  estuary  becomes  thoroughly  mixed  and  any  further  increase  in  w will  have  no 
effect  on  the  density  of  the  fluid  in  the  estuary  and,  therefore,  no  effect  on  qx . 
For  higher  w,  we  have  dqx /&j  = 0.  Therefore,  Eq.  (8. 1)  becomes  independent 
of  u>,  i.  e. , 


■ . = f IJ1  —2  Si  §2 

h ’ h k * h * " * 


(8.2) 


This  tells  us  nothing  about  the  variation  of  qx  with  any  of  the  lengths  of  the  problem 
but  does  indicate  that  the  efflux  is  directly  proportional  to  the  square-root  of  the 
density  difference  between  the  two  fluids  when  the  estuary  is  thoroughly  mixed.  We 
may  extend  this  to  the  more  realistic  case  of  an  outside  linear  density  gradient  and 
find  that  qx  is  proportional  to  the  square  root  of  the  density  gradient. 
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ILet  us  make  a detailed  investigation  of  the  three-layer  circulation  (Long, 

1977).  The  model  is  shown  schematically  in  Fig.  8.  3.  The  estuary- harbor  is 
1 , on  the  left  of  the  figure.  We  assume  that  the  width  of  the  channel  B,  located 

between  the  estuary-harbor  and  the  outside  water,  is  independent  of  z.  The 
, bottom  of  the  channel  and  of  the  two  bodies  of  water  has  an  equation  z = £ (x) 

with  z = 0 chosen  at  a depth  H below  the  level  of  the  outside  body  of  water.  The 
equation  of  the  free  surface  is  z = H + AH  where  AH  is  small  and  only  important 
with  respect  to  the  pressure  distribution.  Outside  of  the  estuary- harbor  there 
is  a two-layer  system  with  an  upper  layer  of  thickness  h.  The  upper  layer  has 

a salinity  S0  and  the  lower  layer  has  a salinity  S2.  The  water  in  the  estuary  is 
thoroughly  mixed  and  has  a salinity  SA.  The  corresponding  densities  are  p0. 

Pi  and  p3.  Later  we  confine  attention  to  a model  with  a uniform  depth  H every- 
where and  a zero  fresh-water  influx.  Here  the  fresh-water  influx  R is  finite 
and  if  the  depth  is  variable  we  require  that  it  be  a minimum  at  the  section  of 
minimum  width  W.  The  fluxes  of  water  of  salinity  S0,  Sj  and  S2  are  q0,  qi 
and  q2,  respectively.  The  lower  interface  has  an  equation  z = z8  and  the  upper 
interface  an  equation  z = z4  so  that  the  thicknesses  of  the  three  layers  are 
(z2~  0 » (zi~zs)  and  (H-ZjJ,  respectively,  if  we  neglect  the  small  quantity  AH 
except  as  it  influences  the  pressure  distribution.  At  the  constriction  the 
thicknesses  of  the  layers  are  T]2H,  (l-lb-ITslH,  and  Tfc,H,  respectively. 

Using  the  hydrostatic  approximation,  the  three  Bernoulli  equations  are 

2 

Pog(H+  AH)  + ■ = fogH  (8.3) 

2B  (H-zx) 

' _2 

(Pi-Po)ezi  + P0S(H+  AH)  + — -SiA — _ - constant  (8.4) 

2B3(Zl-zs)2 


(P2-Pi)gz2  + (Pi-Po)gZi  + Pog(H+AH)  + 


2B3(z3-03 


(P2~Po)g(H-h)  + p0gH 
With  the  Boussinesq  approximation,  these  become 


gAH  + — ao — - = o 
2B  (H-Zl) 


(Ab0-Abi)Zl  + gAH  + “ = constant 

2B2(Zl-z3)a 


z + (Abo-AbJZi  + gAH  + 42 = &>0(H-h)  (8.8) 

2B2(zs-C  )2 

where  we  have  used  the  definitions 

Abj.  ^g^lfij.)  > £b0  =g  ~-p0  (8.9) 

Ps  P2 

Let  us  now  differentiate  Eqs.  (8.  6)-(8.  8)  with  respect  to  x,  evaluating  at  the 
section  of  minimum  depth  and  width  where  9B/3x  = 0,  dc/dx  =0,  B = W,  H-Zl  = TbH, 
zi-zs  = (l-Tb-'ik'.H.Zs-C  = ^2 H.  We  get 


A(gAH)  + —So % = 0 

934  W3T1?H3  9x 


(8. 10) 


~J~(gAH)  + gl-  ■-  

W2  (1— T\j — T)a)3  H3 


^-(gAH)  + (Abo-^bi) 


) = 0 (8.11) 


^(gAH)+  Ab,^+  (Abo-Ab,)^- 


ftza  _ 


w^h3 
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We  may  assume  (Long,  1976)  that  the  interfaces  and  free  surface  do  not  have  a 
zero  slope  at  the  constriction  so  that  the  determinant  of  Eqs.  (8. 10)-(8. 12)  is 
zero.  We  obtain  the  critical  condition 


[(Abo-Abi)-— “ — , j [b1- | - (Abo-AbO  SL 

^(l-n  o-VsfH3^1-  W^fH3-  W2(1-1\3-T12)3H3 

(8.13) 


— “ 1 Abx  - 

W2!}?!!3  L 


W2  T|3  H?  W2(1-T)0-1)2)3I^ 


The  equations  for  conservation  of  volume  and  salt  are 


R + q0  + q2  = qj. 


(8.14) 


qoS0  + qsS2  - qiSj. 


(8.15) 


If  we  use  a linear  relationship  between  density  and  salinity,  we  obtain 


do  = 5 

s r 

= £Lfezlj  R(r-l) 


(8.16) 


(8.17) 


where  s = £b  J &>1,  r = &>0/b2  and  where  ba  is  the  buoyancy  of  water  of  salinity  S2. 
With  the  following  definitions, 

os  a! o3  _ qo2  o _ h 

W2!!3  Ab0  W2H3Ab0  H 


Q|  = — ^ , Q2  = -5 

W2!!3  £b0  W2H3Ab0 


Eqs.  (8.16)  and  (8.17)  become 


Qo  = ^ - ?f  . Qa  = 


s r 


(8.18) 
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Solving  Eq.  (8.23)  for  s,  we  get 

s = m(Qx  , T|o,  T)a)  (8.25) 

Substituting  (8.  25)  into  (8.  24)  we  get 

^3  [Qi » » 71  3 » m(Qi  » 710,  Tla)  1 = P(Qi  * 7)o » T^)  = 0 (8. 26) 

Thus 

Qi  =X(Tlo.Tl3).  s = A(T10, T|a)  (8.27) 

We  assume  that  the  non-dimensional  flux  increases  as  the  mixing  increases. 
When  the  estuary-harbor  is  thoroughly  mixed,  any  further  increase  in  the  mixing 
will  not  change  Qx . As  in  earlier  arguments,  we  assume  that  this  occurs  when 
Qi  is  a maximum,  now  with  respect  to  the  variation  of  two  quantities  71 0 and  tl2 
in  the  equation,  Qx  = x(TI0.  Tfe).  Thus  we  require  for  an  overmixed  estuary 


(8.  26) 


(8.27) 


Mr  ,8-28) 

where  to  make  matters  perfectly  clear,  we  use  subscripts  to  indicate  the  variables 
held  fixed  in  the  partial  differentiation.  To  find  Eqs.  (8.  28)  in  a useful  form,  let 
us  differentiate  (8.  23)  and  (8.  24) 


3*  + 2k 

• T r\m 


+ »k) 

ds  L 


^X.Tla,s  71,  "'(fc.lo.l.  ,!o'71. 

dfx  ax  a^  N a^  . aA  > 

^o.TU.s9^  ^^.Ho.s  98  ^.7,0.^  9^V 


«La) 


ax  ' , af. 

*T\oL 


T1o,*3,s9^o  9Tfe/C 


(8.  29) 


(8.  30) 


+ Sb) 

as  /_  „ _ 

1 Qi  * 71o»  T| a 

9A 

9V-  =0 

(8.31) 

+ *»'\ 

i 3s  /_  _ _ 

8A  , =0 

371a/. 

(8.  32) 

Using  (8.28)  and  eliminating  9A/9T10  and  9A/9T|a,  we  obtain 


Ml 

_ 9fg  9fj 

971 0 

9s 

971 0 9s 

9fi 

9fg 

- Mi 

971 8 

9s 

9 71a  9s 

(8.33) 


(8.  34) 

where  we  now  suppress  the  subscripts.  Eqs.  (8.  33)  and  (8. 34)  together  with  Eqs. 
(8.  21)  and  (8.  22)  are  four  equations  in  four  unknowns,  Qx,  T)0,  Tla,  s(r,  8,  Qf 
given)  and  the  problem  is  determinate.  The  quantities  in  Eqs.  (8.  33)  and  (8. 3'i) 
are 
afi 


= 4Ti0Tlp3s+4in0n2s(s-l)(l-H0)-2Tlon2S(s-l)+2n0[Q1  (s-l)-Qf  f(r-l)  ] 


an. 


(8.  35) 


-4s3  Ho  Til  (1-8) 


f*  =6Tll^s+4TllT)as(s-l)(l-Tlo)'2Tla(Q1-|Qf)a-4s3nlT18(l-P)  (8.36) 

^ = 2Tlltl3+  2713 Tll(s-l)(l-Tl0)+  2T|1 Tils (1  - T]0 )+ 2Tlg  [Qj,  (s-l)-Qf  f (r-1) ] [Qx  - (r-1)  ] 

-n  (8.37) 

+ Qf]fa-4s^lnl(l-8) 

f|Jo  = -3(s-l)a-Tlo-11»)2TllslTl-ls-[Q1(s-l)-Qf |(r-l)Tj  -3(s-l)Q?71§ Tils3 

+ 3 TIq s (S-Dd-Tlo-Tls^^s /lTlls-[Q1(s-l)-Qf|(r-l)]2J+3(Q1-5Qf)2(l-Tl0-Tl8)3Tis 

(8.38) 


-3(1  -710  ■ -D  s)2  (Qi  - f Qf  f T Qt  (s-1) -Qf  f (r-1)  f 
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= ~3(s-l)(l-Tl0-'ns)2T|oS  Tl g s—  TOi  (s-l)-Qf  ^(r-1) ]3  +3T]|Tl|sa  [(s-l)(l-n0-Tla)3-Qf  s] 

-3Tl|TlS(s-l)Q?sa+3(Q1-fQf)2(l-'n0-na)8Tl38-3nls(l-n0-n3)3(Q1-|Qf)a  (8.39) 

-3(1-Tl0-T)a)s (Qx-f  Qf)2[Qi(s-l)-Qf  f (r-l)]3+  SQ^Tl^Q,-  | Qt)a 

= (l-Tlo-'n2)3-Qf  ^Is-CQ^s-D-Q,  f(r-l)]2  Tlgs-(4s3-3s2)Qfn037ll 
+ T,o  (s-Dd-no-nZ-Qfs  nls-CQj  (s-l)-Qf  | (r-1) ]2  , 

+ n3s  (s-1)(1-T10-T18)3-Q2s  nl-ZCQ^s-D-Q,  ffr-DJCOi-Q, 

(8.  40) 

+ 2^(Qi-  |Qf)nls(l-n0-Tls)3-  Tl38(l-Tlo-Tl8)3(Ql-f  Q'f 

-2  (Qi - f Qf ) [Qi  (s-l)-Qf  | (r-1) ]S (1 — T| o — T| 2 )3  -2QX (Qx - | Qf)  ^fs3Tfs 

+ 2 [Qj.  (s-l)-Qf  f (r-1) ] [Qi  - ^-(r-1) ] [Qx  - § Qf  f (1-T)0  -H s)3  + 2Q3  si)3  (Qx  - f Qr  f 

The  solution  of  the  general  problem  involves  the  solution  of  the  four  complicated 
simultaneous  algebraic  equations  and  we  will  limit  attention  to  the  special  case  in 
the  next  section.  We  note,  however,  that  for  a zero  fresh-water  influx,  the  solution 
has  the  simple  form  Qx  = f(0)  so  that  the  discharge  qx  is  proportional  to  the  square 
root  of  the  density  difference  between  the  two  outside  fluids  as  we  found  out  from 
dimensional  analysis. 

Let  us  now  consider  the  special  case  in  which  the  whole  system  has  a uniform 
depth  H,  the  fresh-water  flux  is  zero,  and  the  two  fluids  outside  of  the  harbor  have 
equal  thicknesses.  Then  0 = i and  from  symmetry  s = 2 and  = 710.  Eq.  (8.  21) 
becomes  an  identity  and  Eq.  (8.  22)  becomes 
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4T|o(1-2T10)3-  l6QtTl06-  4 n 3(1  - Z 71  <^3  + 8Qj T1 3 + Qt (1  - 2 TR0 )3  = 0 (8.41) 

The  condition  for  overmixing  is  90^/ 971 0 =0.  Differentiating  Eq.  (8. 41)  we  get 


-12n®(l-2Tl0)s  + 12Q127)03(l-2'n0)a-3Q?(l-2-Ib)2+12Tf0(l-2Tl0): 


(8.  42) 


-48Q^  710s  + 12Q14n^-6QfTl2(l-2H0)3  =0 


The  two  equations  (8. 41)  and  (8. 42)  determine  the  problem  and  the  relevant  solu- 


tion is 


"Ho  =4.  of  = 


(8.43) 


In  dimensional  terms,  the  outgoing  fluid  has  a thickness  of  H/2  and  a flux 


XX  »»  IX. 

qi  =— j £ 


(8.  44) 


A laboratory  model  was  constructed  to  study  a three-layer  circulation.  A 
long  channel  of  width  W0  was  used  as  shown  in  cross-section  and  in  plan  view 
in  Fig.  8. 4.  The  opening  between  the  harbor  and  the  outside  region  has  a uniform 
width  W.  The  outflowing  fluid  has  a vertical  thickness  aH  far  downstream.  Atten- 
tion was  confined  to  the  case  of  fluids  of  equal  depth  so  that  0 = J. 

Let  us  first  consider  the  steady-state  theory  for  an  infinitely  long  channel. 

A portion  of  the  problem  is  identical  to  that  considered  above  and  we  are  led  to 
equations  identical  to  Eqs.  (8. 41)  and  (8. 42) . The  solution  is  again  7)0  = | , 

Qi  = 1/64.  To  obtain  an  equation  for  a,  we  use  the  Bernoulli  equations  for  the 
outflowing  fluid  and  for  either  the  upper  or  lower  fluid.  These  are 


gAH  + 


2W2  71o  H2  W^H^l-a)" 


-f 0 (l-Tlo)H  + gAH  + ^ = H(i+a)  + 

2W3  (1— 2Tlo)3  H3  4 2W;aaHa 


(8. 45) 


(8. 46) 
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Eliminating  gAH  between  these  two  equations,  we  get 

1 Up  a Q?  w3  Q3  Q?  _ Cjf  W2 

4 ' 2 " 4 + 2(l-a)a  V^o  ‘ 8%  + 2(1-2T>d)2  ' 2a2  W<, 


= 0 


(8.  47) 


Using  Tq  = l/4,  Q2  = l/64  and  imposing  the  condition  a = 0,  W/W0  s=  0,  we  obtain 

a = |[1-(1- (8.48) 

A graph  of  this  relationship  is  shown  in  Fig.  8. 5. 

The  experiment  was  conducted  by  putting  a two-layer  system  of  equal  depths 
in  the  channel  and  then  using  several  egg-beaters  to  mix  the  fluid  in  the  harbor.  The 
fluid  in  the  harbor  was  thoroughly  mixed  in  a few  seconds  and  the  middle  layer  began 
to  move  down  the  channel  in  a surge.  Evidently  a quasi-steady  state  resulted  in  the 
vicinity  of  the  mouth  of  the  harbor  before  the  surge  was  able  to  reach  the  other  end 
of  the  channel  and  reflect.  The  middle  layer  was  unfortunately  still  quite  turbulent 
near  the  mouth  and  for  some  distance  along  the  channel,  but  to  the  extent  that 
observation  of  its  thickness  was  possible,  it  was  close  to  H/2  at  the  mouth  in 
accordance  with  the  theory  and  as  evidenced  by  the  photograph  of  Fig.  8.  6.  The 
thickness  ratio  a was  considerably  larger  than  Eq.  (8. 48)  indicates,  probably 
because  of  the  presence  of  a "drowned"  hydraulic  jump  downstream. 
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